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Summary 

The focus of this research is to numerically predict an 
infrared image of a jet engine exhaust plume, given field 
variables such as temperature, pressure, and exhaust 
plume constituents as a function of spatial position within 
the plume, and to compare this predicted image directly 
with measured data. This work is motivated by the need 
to validate Computational Fluid Dynamic (CFD) codes 
through infrared imaging. The technique ot reducing 
the three-dimensional field variable domain to a two- 
dimensional infrared image invokes the use of an inverse 
Monte Carlo ray trace algorithm and an infrared band 
model for exhaust gases. 

This report describes an experiment in which the above- 
mentioned field variables were carefully measured. 

Results from this experiment, namely tables o( measured 
temperature and pressure data, as well as measured inlra- 
red images, are given. The inverse Monte Carlo ray trace 
technique is described. Finally, experimentally obtained 
infrared images are directly compared to infrared images 
predicted from the measured Held variables. 

1. Introduction 

1.1 Motivation 

The validation of Computational Fluid Dynamics (CFD) 
codes is important in the design of future generations of 
both military and civilian aircraft. Infrared imaging can be 
used as a validation tool for CFD codes, particularly for 
Short Take-Off and Vertical Landing (STOVL) aircraft. 
The validation technique involves a quantitative com- 
parison between an experimental inlrared image and a 
predicted infrared image, with the latter based on the CFD 
solution. This report describes the assembly of a data base 
which has been completed in support ot an effort to 
validate CFD codes through infrared imaging. An experi- 
ment has been conducted in which the flow field of a 
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modified auxiliary power unit (APU) has been surveyed, 
where temperature and pressure data, as well as corre- 
sponding infrared images have been obtained. Infrared 
images predicted from this data base arc compared to 
experimentally obtained infrared images, thereby 
demonstrating an important step in the CFD code 
validation technique. 

In recent years improved algorithms and advanced high- 
speed computers have led to increasingly complex CFD 
analyses. Complete Navier-Stokes analysis of such 
complex aerospace vehicles as the National Aeronautics 
and Space Administration’s (NASA) Space Shuttle 
(ref. 1 ), the Stratospheric Observatory For Infrared 
Astronomy (SOFIA) Airborne Observatory (ref. 2), and 
the McDonnell Douglas Harder YAV-8B (ref. 3) have 
been successfully performed. However, validation of the 
predicted results for these complex acropropulsive flow 
fields has been hampered by the lack of reliable test data. 
For example, definition of the complex jet-induced 
interactions produced by the Harrier operating in ground 
effect, obtained using conventional measurement 
techniques such as pressure transducers and thermo- 
couples, has proven impractical due to the si/e scale and 
dynamic nature of the flow field. In addition, a critical 
need exists for validated CFD tools to ensure the 
successful design of nozzles and propulsion systems lor 
future aircraft. 

The motivation for the interest in STOVL aircraft stems 
from a requirement to operate aircraft, especially fighter 
aircraft, from a minimum operating airstrip or bomb- 
damaged runway. Interest in short-field takeoff capability 
extends from large land-based aircraft to aircraft carrier 
launch operations and ground-based remote forward-area 
operations. Current STOVL aircraft include vehicles such 
as the McDonnell Douglas Harrier AV-8B (ref. 4). 

The aerodynamic flow field about a Harrier AV-8B 
operating in ground effect is difficult to model. In this 
flight regime, lift is augmented by four rotatable nozzles 
positioned fore and aft of the aircraft’s center of gravity. 
Also, small reaction jets arc located at the aircraft’s nose, 
tail and wing tips. The ground effect flight regime is 



characterized by a small forward velocity component and 
regions of high subsonic to supersonic jet flows. The 
recirculating fluid dynamic interaction between the 
vectored jets, the ground, and the airframe can cause 
several effects including propulsion-induced loss of lift, 
or “suck down,” ingestion of foreign objects or debris 
(FOD), and ingestion of warm, oxygen depleted, or 
vitiated, air (refs. 5-7). Attempts to model these effects 
by conventional means, such as wind-tunnel experiments, 
have been unsuccessful, with the end result being that 
full-scale powered tests are required to study the actual 
flow Held. 

The United State Air Force would like to replace the 
aging F- 1 6 aircraft fleet. Also, the Navy would like to 
replace the A-6 and the Marines would like to replace 
the Harder AV-8B. In light of shrinking defense budgets, 
the United States government is considering a multirole 
tactical airplane to meet the requirements of the Air 
Force, Navy, and Marines (ref. 8). The US Advanced 
Research Project Agency (ARPA) hopes to award a 
STOVL Strike Fighter (SSF) contract in the late 1990s. 
Major airframe and engine companies are working now in 
hopes of making a successful bid for the ARPA contract. 
For example, wind-tunnel testing of an SSF model with 
lifting fan, made by British Aerospace and McDonnell 
Douglas, is planned for 1995 at NASA Ames (ref. 4). 
NASA Ames has been working in support of this effort in 
a wide variety of areas, particularly in the area of STOVL 
technology (ref. 9). Research activities include wind- 
tunnel testing of small-scale models (refs. 10-12), CFD 
modeling of critical flight regimes (refs. 13 and 14), and 
control and stability analysis in the Vertical Motion 
Simulator (VMS) (refs. 15 and 16). 

The Powered-Lift Technology Branch at NASA Ames 
is directly involved in a number of STOVL initiatives. 

In support of their work in this area, an Agema 
Thermovision® Dual 880 Infrared Imaging System 
has been purchased in an effort to investigate the 
feasibility of using infrared imaging to validate CFD 
codes. The imaging system employs two scanning 
infrared cameras that are sensitive in the medium 
wavelength band (MWB), 2.5 to 5.5 Jim, and in the long 
wavelength band (LWB), 8 to 12 |im. The infrared 
detectors arc cryogenically cooled and provide a source 
temperature sensitivity on the order of 0. 1 °C at 30.0 °C. 
The system operates by scanning a scene and digitally 
recording the infrared image to a hard disk (ref. 17). A 
thorough technical discussion of the Agema infrared 
imaging system is available from the Agema company 
(ref. 18). Also, the imaging system is discussed further in 
Section 2 of this report. 


Since its delivery to NASA Ames in November of 1989, 
the Agema infrared imaging system has been useful as a 
nonintrusivc measurement technique which provides real- 
time flow Field images of sufficient detail to resolve the 
small- and large-scale flow structures and surface heating 
effects needed for the applications described above. The 
practicality of using infrared imaging as a means of 
examining these types of complex aeropropulsive flow 
fields has been demonstrated in a series of tests conducted 
at NASA Ames. Birckelbaw and Nelson (ref. 19) have 
shown the ability of an Agema system to resolve small- 
and large-scale flow structures, airframe surface heating, 
and interactions between the ground and plume. Inflight 
and ground-based infrared images of the Harrier YAV8B, 
the XV- 1 5 Tiltrotor Aircraft, and the E-7 STOVL 
research model were shown to have value for the 
qualitative examination of the flow fields. 

An example of the type of results which can be obtained 
from the infrared imaging system is shown in figure 1 . 
Figure 1(a) shows an inflight photograph of the 
McDonnell Douglas Harrier YAV-8B and figure 1(b) 
shows a MWB experimental infrared image of the Harrier 
in a 30-ft (9.14-m) altitude and 30-kt ( 15.4-m/s) velocity 
flight condition. Because of the strong emission in the 
infrared, the hot-exhaust stream surrounding the aircraft is 
clearly visible. Also, the effects of ground heating and 
aircraft surface heating are visible. The “hot spots” which 
are visible on the fore and aft of the aircraft are exhaust 
ports, sometimes called “puffers.” These “puffers” are 
used to maintain aircraft stability. 

While the usefulness of infrared imaging as a How- 
visualization tool is being explored, there exists a 
concurrent effort to better understand infrared imaging 
and modeling of exhaust plumes and gaseous radiation. 
Over the past six years joint research at NASA Ames and 
the Virginia Polytechnic Institute and State University 
(VPI&SU) has been underway to develop a compu- 
tational technique which would permit quantitative use of 
measured infrared data for validation of advanced CFD 
codes (refs. 20-24). This research has concentrated on the 
development and validation of an inverse Monte Carlo 
radiative analysis technique that accurately generates a 
two-dimensional infrared image directly from an existing 
three-dimensional CFD solution. 

The inverse Monte Carlo ray-trace algorithm was recently 
completed by Turk and co-workers (ref. 25). Figure 2 
illustrates how the algorithm works. The computer code 
traces a ray from the detector of an infrared camera, 
through the CFD solution space, to some destination cell. 
The location of the destination cell is dependent upon 
the initial orientation of the ray, the absorptance of the 
destination cell, and the absorptance along the path from 
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(he camera objective to the destination cell. The radiative 
contribution of each ray is computed based on the 
properties of the path traversed by the ray, and by the 
properties of the destination cell. Creating the two- 
dimensional infrared image involves tracing a large 
number of these rays for each pixel on a virtual two- 
dimensional image screen. Increasing the number of 
pixels and increasing the number of rays per pixel 
increases the accuracy of the predicted image, within the 
limits permitted by the resolution the CFD solution. 

Development of this technique began with a proof-of- 
concept experiment by Hardman (refs. 20 and 21). In 
Hardman’s experiment, an infrared image of an exhaust 
plume issuing from an APU was measured and compared 
with an image computed from the experimentally 
measured temperature distribution. A gas turbine-driven 
APU and an attached flow straightening assembly were 
used to create a uniform heated jet. A stepper-motor- 
driven combination Kiel and thermocouple probe was 
used to measure the pressure and temperature profiles at 
regular intervals downstream of the nozzle exit. Although 
the thermocouple probe was not corrected for conduction 
and an overly simple water-vapor-only band model was 
used to convert the measured temperature field into an 
infrared image, the experimental and predicted infrared 
images compared favorably. The results demonstrated the 
concept that infrared images could be generated from 
knowledge of the field variables and compared with 
experimentally obtained infrared images. 

Hardman’s experiment and results were successful and 
important enough to foster a continued research program 
in this area. Also, it was of interest to extend the ray-tracc 
technique to the area of CFD code validation, particularly 
for STOVL applications. The next logical step to real- 
izing the goal of predicting infrared images which can be 
compared quantitatively to experimental infrared images 
was to improve the water-vapor-only band model used in 
Hardman’s computer program. 

Over the past twenty years, considerable research has 
been done on infrared band models (refs. 26-30). Several 
commercial codes refs. 3 ! and 32) and a few company 
proprietary codes (ref. 33) are available. Each band model 
has its particular advantages; for example, some are best 
for high-temperature applications. In concert with the 
research objectives described above, Nelson (ref. 22) 
describes the selection of a public-domain infrared 
gaseous radiation band model which has been tailored for 
CFD code validation. In reference 22, Nelson describes 
the NASA Band Model, the band model published in 
NASA SP-3080, the Handbook of Infrared Radiation 
from Combustion Gases (ref. 34), and its improvement by 
implementing data tables from the North Atlantic Treaty 


Organization (NATO) Infrared Air Target Model 
(NIRATAM) (ref. 32). The improved band model 
includes the effects of water vapor, carbon dioxide, and 
carbon monoxide, the three gases which emit most 
strongly in the infrared. The improved band model 
provides constituents at a five-wavenumber resolution, a 
sufficiently high resolution for most CFD code validation 
applications. 

The improved infrared band model has been incorporated 
into an inverse Monte Carlo ray-tracc algorithm. The 
algorithm is similar to Hardman’s code in that it reduces a 
three-dimensional solution field to a two-dimensional 
infrared image. However the latest algorithm, recently 
completed by Turk (ref. 25), operates directly on a 
generalized overset grid topology and CFD solution set. 
This means that the Monte Carlo technique, as it has been 
coded in the most recent computer program, can be 
applied to a variety of CFD solution sets. 

A qualitative visual comparison between figures 3 and 4 
indicates the most recent results from Turk's inverse 
Monte Carlo ray-trace code. Figure 3(a) shows a 
photograph of the Boeing 747 Space Shuttle Carrier 
Aircraft and Figure 3(b) shows an experimental inflight 
infrared image of the same Boeing 747. The two red and 
white exhaust plumes, shown in figure 3(b), correspond 
with the inboard and outboard engines of the aircraft. 

Also, the reflected energy can be seen on the underside of 
the wing in figure 3(b). Notice that figure 3(b) is a MWB 
infrared image; the image is taken over the band of 2.5 to 
5.5 (im. 

Figure 4 is an infrared image which has been predicted 
from an existing CFD solution set about the Boeing 747 
in a similar flight condition. A qualitative comparison 
between figure 3(b) and figure 4 reveals that the Monte 
Carlo algorithm is useful in identifying objects such as 
the fuselage, and for locating the size and structure of the 
exhaust plumes. The predicted image shown in figure 4 is 
a spectral image; the image is predicted at 4.3 pm. While 
a quantitative comparison between a broad-band image, 
such as the image shown in figure 3(b), and a spectral 
image, such as the image shown in figure 4, is not 
possible, the qualitative comparison between the two can 
be useful, particularly in identifying areas which may 
need improvement in the predicted CFD solution. The 
Monte Carlo algorithm and issues pertaining to infrared 
image prediction are discussed further in Section 4 and in 
reference 25. 

1.2 Goal 

Before the Monte Carlo algorithm can be used with 
confidence to validate a CFD solution, it too must be 
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validated. Such a process begins by generating a 
relatively simple CFD-like solution from a known 
temperature and pressure field. Next the Monte Carlo 
technique is used to predict an infrared image from the 
CFD- like solution. Finally, the predicted infrared image 
is compared with an experimentally obtained infrared 
image. It is the goal of the research described in this 
report to assemble the data base required for such a 
comparison. 

Creating the CFD-like solution involves overlaying a 
known temperature and pressure field onto a CFD-like 
grid. The CFD-like grid is necessary because the Monte 
Carlo technique has been coded to directly interrogate a 
CFD solution. The data set that serves as the source for 
the temperature and pressure measurements must also 
include corresponding infrared images. The required data 
base of known temperature and pressure distributions, 
with corresponding infrared images, is not generally 
available. 

This report describes an experiment which has been 
conducted in support of the research efforts described 
above. The experiment involves repeating Hardman’s 
(refs. 20 and 21 ) experiment by conducting a complete 
temperature and pressure survey of the How field exiting 
a gas turbine-driven AFU. The experiment is unique and 
different from Hardman’s work in that temperature and 
pressure surveys were taken simultaneously with corre- 
sponding infrared images. Also, efforts were made to 
ensure the accuracy and repeatability of the flow field 
surveys. The experiment was conducted in an effort to 
provide a data base from which the inverse Monte Carlo 
ray-trace technique and a potential CFD solution of this 
flow could be validated. Infrared images, temperature and 
pressure plots, velocity plots, and meteorological data 
from the experiment arc given. Finally, infrared images 
predicted from these data using the inverse Monte Carlo 
ray-trace code arc presented and compared to experi- 
mentally obtained infrared images. 

2. Experiment 

An experiment has been conducted in which an exhaust 
plume issuing from a modified auxiliary power unit 
(APU) has been surveyed. Exhaust flow exiting from both 
an axisymmetric and a rectangular nozzle were surveyed; 
however, only the axisymmetric data are given in this 
report. Temperature and pressure data were obtained 
using standard temperature and pressure probe tech- 
niques. Infrared images were obtained using a scanning 
infrared camera. The APU experiment is described in this 
section. 


2.1 Background 

In June 1989, Robert Hardman, a graduate student in 
Mechanical Engineering at VPI&SU, began work at 
NASA Ames with funding from NASA’s postbacca- 
laureate program. One of Hardman’s objectives during 
his eleven-month stay at NASA Ames was to prove the 
concept that infrared images could be predicted from a 
known temperature and pressure field. As a result of his 
work, which is thoroughly described in his Master of 
Science thesis (ref. 20), an APU was converted to provide 
an axisymmetric heated exhaust flow. Specifically, 
Hardman designed and attached a stainless steel flow 
straightening assembly and converging nozzle to the 
exhaust port of an APU (described in Sections 2.3.5 
and 2.3.6). Hardman surveyed this exhaust flow with a 
combination thermocouple and Kiel probe. The results 
from this temperature and pressure survey were success- 
fully used to predict an infrared image. Hardman 
completed his Master’s degree in May 1990. 

One month after Hardman graduated, Ed Nelson started 
graduate school in Mechanical Engineering at VPI&SU. 
During the following summer, in June 1991, Nelson 
started work at NASA Ames under the same post- 
baccalaureate program as Hardman. The focus of 
Nelson’s work was the development of an infrared band- 
model which could be used for computational fluid 
dynamic (CFD) code validation applications; this work 
is described in his Master of Science thesis (ref. 22). 
Nelson’s work was motivated by the need to improve 
upon the “weak link” in Hardman’s infrared image 
generation code, the infrared band model. 

Before Nelson returned to VPI&SU in January 1992 to 
complete his Master’s degree, discussions between 
VPI&SU and NASA Ames revealed a need to develop a 
computer program which could predict infrared images 
directly from CFD solutions. Such a computer program 
would have applications for CFD code validation as well 
as infrared signature prediction for Department of 
Defense (DOD) applications. The computer program 
would be similar in concept to Hardman’s code, but 
written differently to provide the user with options of 
image quality and image accuracy. An agreement was 
reached between Dr. J. Robert Mahan, a professor of 
Mechanical Engineering at VPI&SU, and Dr. Larry 
Birckelbaw, an aerospace engineer and civil servant at 
NASA Ames, that a graduate student from VPI&SU 
would begin work on this project. 

In June 1992, Jeffrey Turk, a graduate student in 
Aerospace Engineering at VPI&SU, began work at 
NASA Ames. Turk’s work focused on the development 
of a computer program which could predict infrared 
images from CFD solutions. The resulting computer code 
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incorporates a Monte Carlo-based ray-trace algorithm 
which can predict an infrared image directly from a CFD 
solution. This computer program is described further in 
Section 4 and in reference 25. During the first six months 
while Turk worked to develop this computer program, the 
need for a data base to validate this algorithm became 
apparent. 

In January 1993, Nelson returned to NASA Ames with 
the objective of creating a data base which could be used 
to validate CFD codes, as well as to validate Turk’s 
program for predicting infrared images. For these 
applications it was necessary that such a data base include 
a detailed temperature and pressure mapping of a heated 
exhaust flow. Corresponding infrared images of the 
exhaust flow were also necessary to complete the data 
set. By virtue of the accessibility to an infrared imaging 
system and the APU previously modified by Hardman, 
NASA Ames became the logical choice for an experiment 
site. 

During the following fourteen months, from January 1993 
to March 1994, Nelson worked to complete this data base. 
After acquiring the necessary experimental equipment 
and the personnel trained to operate that equipment, a 
preliminary experiment was conducted. On Thursday and 
Friday, May 14 and 15, 1993, Dr. Mahan visited NASA 
Ames to participate in the preliminary experiment. Six 
experimenters, including Dr. Mahan and Dr. Birckelbaw, 
tested the techniques required to conduct such an 
experiment. 

From the preliminary experiment in May 1993, we 
learned that the outdoor testing should be conducted at 
night to avoid the high winds that are common in the 
summer afternoons of the San Francisco Bay Area. 
Concerning the acquisition of infrared images, testing at 
night offered the advantage of a decrease in background 
radiation over daytime testing. Also, we recognized the 
importance of leaving our equipment set up for the test; 
the assembly of the equipment required approximately 
sixteen man-hours of labor. Electrical interference 
between the power supply, bugs in the program tor data 
acquisition, the requirement for a better backdrop tor 
infrared data, the convenience of measuring atmospheric 
conditions, and personnel comfort were all factors that 
needed improvement for future tests. Several post- 
experiment meetings involving discussions among the 
experimenters on how to improve the process resulted in 
a much-improved experimental procedure. 

During the week of July 12-16, 1 993, the APU exhaust 
plume issuing from an axisymmetric nozzle was 
surveyed. The axisymmetric nozzle is shown as part of 
the flow straightening assembly which is described in 
Section 2.3.6. During the eight work days between 


October 18 and 27, 1993, the APU exhaust plume issuing 
from a rectangular nozzle was surveyed. Data from both 
APU experiments will be published in NASA Technical 
Memorandums (NASA TMs). 

The rectangular nozzle data arc not described in this 
report because of space limitations. One goal of this 
report is to show how the measured temperature and 
pressure data can be reduced to a CFD-likc data file and 
subsequently used to predict an infrared image which is 
comparable to an experimentally obtained infrared image. 
That comparison is available only for the axisymmetric 
data set at this time. While the data from the rectangular 
nozzle experiment has been reduced, it has not yet been 
converted into a CFD-like data file. Once the CFD-like 
file has been created, some lead-time is necessary to de- 
bug the infrared image prediction computer program used 
to analyze the data file, a potentially time-consuming 
task. For these reasons only the axisymmetric data are 
presented here. 

2.2 Facility 

The experiment was conducted at the Outdoor Aero- 
dynamic Research Facility (OARF) of the National Full- 
Scale Aerodynamics Complex (NFAC) at NASA’s Ames 
Research Center at Moffett Field, California. A trailer and 
underground control room provided excellent storage for 
equipment such as computers and the data acquisition 
system. Also, these rooms provided personnel safety and 
minimum structure interference during the outdoor test. A 
1 00- ft (30.5-m) square concrete pad provided an excellent 
place to conduct the experiment. A weather station was 
also located at the facility. Jet engine fuel, electricity, 
lighting, and safety equipment such as fire extinguishers 
and 91 1 service were also available at the OARF (ret. 35). 

2.3 Description of Test Equipment 

Figure 5 is a photograph taken during the APU experi- 
ment at the OARF. Equipment such as the aluminized 
Mylar backdrop, the large and small traverses, the APU, 
and combination temperature and pressure probe can be 
seen in the photograph. These items and other test 
equipment arc described below. 

2.3.1 Large Two-Axis Traverse- A large two- 
dimensional traverse, shown in figure 6, was used to 
position the probe and to support a smaller three- 
dimensional traverse rig. The traverse provided a maxi- 
mum horizontal displacement of 40 in. ( 101 .6 cm) and a 
maximum vertical displacement of 40 in. (101 .6 cm). 
Covering an area of approximately 30 ft 2 (9. 14 m 2 ) and 
having a mass of approximately 1 200 lb (544 kg), the 
large traverse provided a very stable support structure lor 
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other test equipment such as the smaller traverse and 
probe. A few modifications were made to the existing 
traverse rig specifically for this experiment. Namely, all 
of the wooden components were replaced by aluminum 
components, an improvement that was necessary to meet 
fire code safety requirements. The vertical and horizontal 
position adjustments to the large traverse were made 
using a Shaw Engineering motor-controller which is 
described in Section 2.3.4. 

2.3.2 Small Three-Axis Traverse- A smaller three- 
dimensional traverse was “bolted” to the large traverse 
described in the preceding section. The smaller traverse 
was used to position the combination temperature and 
pressure probe within the exhaust stream. This traverse 
provided a freedom-of-movement of 18 in. (45.7 cm) 
in all three directions. The traverse could position the 
probe within the flow with an accuracy of ±0.001 in. 
(±0.025 mm) in all three directions. A Velmex motor- 
controller was used to operate the three stepper-motors of 
the traverse. Figure 7 is a photograph of the small three- 
dimensional traverse. Here, the small traverse is shown 
attached to the larger traverse for the experiment. The 
small traverse rig and the corresponding motor controller 
were on loan for this experiment from the Naval Post- 
Graduate School in Monterey, California (ref. 36). 

2.3.3 Motor Controller for the Small Traverse- 

A Velmex 8300 Series Stepping Motor Controller/Driver 
was used to operate the small three-dimensional traverse 
(ref. 37). The controller utilized a RS-232C connector that 
accepted ASCII characters as input. The baud rate was 
switch-selectable between 50 and 9600 BPS. The Velmex 
controller was connected to the stepper motors of the 
small traverse by three relatively short cables. Because of 
the short cables, the controller was positioned on the shelf 
of the large two-dimensional traverse alongside the three- 
dimensional traverse. A digital display on the face plate 
of the controller showed the distance traveled by the 
probe from a pre-set home position. The motor controller 
could be operated directly from the face plate, or remotely 
from the trailer via a computer and cable. 

2.3.4 Motor Controller for the Large Traverse- 

A Shaw Engineering motor-controller was used to posi- 
tion the large two-dimensional traverse rig. The controller 
could position the probe with an accuracy of ±0.001 in. 
(±0.025 mm) (ref. 38). A digital display on the face plate 
of the controller showed the distance traveled from a pre- 
set home position. Unlike the controller for the small 
traverse rig, this controller could not be operated remotely 
from computer. Adjustments to the large traverse were 
made at the controller itself, requiring someone to leave 
the test trailer and make traverse adjustments in a noisy 
test environment. Also, the motor-controller was designed 


to send output to only one motor. Consequently, two such 
controllers were required to operate the two-dimensional 
traverse rig. However, only one operational controller was 
available during APU experiment. If the large traverse 
needed to be moved in two directions, the cables leading 
from the controller to the motor had to be switched. 

2.3.5 Modified Auxiliary Power Unit— An FAA 

Category II, Class B Gas Turbine Auxiliary Power Unit 
(APU) was utilized to produce the hot free -jet, or plume. 
The APU is an Airesearch Model GTCG30- 142C 
manufactured by Garrett Auxiliary Power Division of the 
Allied Signal Aerospace Company (ref. 39). Attached to 
the nozzle exit is a flow-turning pipe, a diverging nozzle, 
and a flow-straightening assembly, which are described 
below. The APU was originally designed as an air-starter 
for larger jet engines. Modifications to the APU, namely 
the addition of a ten-gallon fuel tank and the flow 
straightening assembly, were made by Hardman (refs. 20 
and 21). Figure 8 is a photograph of the APU in which the 
nozzle assembly and fuel tank can be seen. 

The APU requires a continuous source of 28-Vdc power. 
This power was supplied by a 480-Vac-to-28-Vdc 
rectifier borrowed from the NFAC. Figure 9 is a photo- 
graph of the Hobart rectifier (ref. 40). Use of the rectifier, 
as opposed to other power-supply options such as a 
portable large-aircraft air-starter-system, decreased noise 
at the test site and reduced electrical interference to the 
data acquisition system. 

The APU used JP-4, a jet engine fuel used in practically 
all European and U.S. military jets (ref. 41 ). The ten- 
gallon fuel tank was replaced, for this experiment, by a 
direct fuel line to a large 10 ,000 gal (37,854 I) lank at the 
OARF. Tapping into the large tank permitted the APU to 
be operated continuously, an improvement over earlier 
tests which offered a substantial time savings during the 
APU test. The APU start- and kill-engine switch was 
located on an attached instrument panel, along with an 
exhaust gas temperature gauge. An occasional visual 
inspection of the temperature gauge was required to make 
sure the APU did not overheat during continuous 
operation. 

2.3.6 Flow Straightener- A flow-straightening section 
and axisymmetric exhaust nozzle assembly were added 
to the APU by Hardman (ref. 20). The straightening- 
assembly, shown in figure 10, includes a clamping flange 
and 90-deg elbow, a Zanker tube bank, an eddy-removu! 
screen, and an axisymmetric nozzle. The assembly is 
constructed entirely from 20 gauge (0.035 in.) 304 
stainless steel. An adjustable flow collar, designed to 
alleviate the additional back-pressure induced by the 
assembly, is installed immediately downstream of the 
90-deg elbow. No problems, such as engine stalls due 
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to the added back-pressure from the assembly were 
encountered. For this reason the flow collar remained 
completely closed for all APU tests. 

The Zanker tube bank, shown in Figure 1 1, is positioned 
downstream of the 90-deg elbow and adjustable flow 
collar. The tube bank helps to straighten the exhaust 
flow (ref. 42). An eddy-removal screen positioned 
downstream of the tube bank serves to eliminate large- 
scale eddies. The stainless steel mesh has 450 wires per 
inch (177 wires/cm) in each direction of the screen. The 
design and fabrication of the straightener is discussed 
fully by Hardman (ref. 20). 

2.3.7 Combination Kiel and Thermocouple Probe- 

A combination Kiel and thermocouple probe, shown in 
figure 12, was used to measure the total pressure and total 
temperature of the exhaust stream exiting the APU. The 
probe is 1 2.8 in. (32.5 cm) long and has a diameter of 
0. 1 88 in. (0.478 cm). The length and diameter of the 
probe were selected to provide adequate stiffness. 

Because of drag forces, stiffness was an important 
concern when attempting to maintain the probe's position 
in the flow; it was desired to keep the jet from deflecting 
the probe. A large aluminum brace was added to the base 
of the probe to increase its stiffness for this reason. 

The Kiel probe was manufactured by the United Sensor 
Division of the United Electric Controls Co.; this 
company has since changed its name to the United 
Electric Controls Co., Flow Sensor Products Division 
(ref. 43). The yaw and pitch angles, defined in figure 13, 
of the tlow entering the type-KA Kiel probe could vary 
up to ±50 deg without error in the total pressure reading. 

A true total pressure is indicated up to a Mach number of 
1.0. Also, the probe is insensitive to Reynolds number 
except at extremely low velocities; the limiting velocity is 
approximately 4.0 ft/s. Errors introduced by turbulence 
are negligible, except for severe turbulence which may 
decrease slightly the acceptable range oflhe yaw and 
pitch angles (ref. 43). The typc-KA Kiel probe was 
selected for the range of velocities encountered during the 
APU experiment. It is estimated that the yaw and pitch 
angles did not exceed 5 deg over the duration of the test. 

The type-K, Chromel-Alumel, thermocouple is located 
0. 155 in. (3.96 mm) from the Kiel probe head. A second 
typc-K thermocouple was used to monitor the ambient air 
temperature. Both thermocouples were connected to a 
signal conditioning unit, discussed in Section 2.3.10, 
which housed the thermocouple modules. 

2.3.8 Thermocouple Calibration Unit- Calibration of 
the thermocouples was conducted using an Omega Hand- 
Held Calibration Unit, Model CL23-JKT Calibrator 
(ref. 44). The Calibration Unit checked the accuracy of 


the thermocouple by comparing the output voltage from 
the thermocouple to a reference voltage. The Calibration 
Unit provided an accuracy of 0. 1 percent of the recorded 
temperature plus I °F. For example, if the recorded 
temperature was 800 °F, that measurement would be 
accurate to 1 .8 °F. 

2.3.9 Kiel Probe Calibration Unit- A pressure 
calibration unit (PCU) was used periodically during the 
experiment to test the response of the Kiel probe. The 
PCU was a DPI-605 manufactured by Druck Inc. 

(ref. 45). The unit was quoted as having an accuracy of 
0.05 percent of full scale (30 psid max.), which implies 
accuracies of 0.015 psid. The unit was attached to the 
pressure lines with Swage fittings. A pressure-relief valve 
on the PCU was opened to “zero” the pressure reading. 
Next, the relief valve was closed and a pressure was 
applied. The pressures measured by the PCU and the 
pressure measured by the computer system were 
recorded. After a number of points, these values were 
compared. The process was repeated for each transducer. 
The PCU itself was calibrated just one week prior to the 
test. 

2.3.10 Signal Conditioning Unit- A signal conditioning 
box (SCU), shown in figure 14, was constructed to house 
a series of analog signal conditioning modules. In the 
photograph of figure 14, the lop cover of the SCU has 
been removed to show the linear and nonlinear 
thermocouple modules. These modules. Models 5B47, 
were manufactured by the National Instruments Co. 

(ref. 45). The modules have excellent noise resistance 
and amplification for millivolt sources. For the APU 
experiment, only the linear modules were used; these 
modules provided consistent treatment of the signal 
arriving from the thermocouples, unlike the nonlinear 
modules which were tried and gave erroneous and 
inconsistent results. The low-level input signal, in the 
millivolt range, arriving from the thermocouples was 
amplified to between 0 and +5 Vdc. The amplified 
voltage directly corresponds to a temperature between 

0 and 500 °C. The modules were fully encapsulated and 
required no field adjustments. 

2.3.11 Data Acquisition System- A Macintosh computer 
equipped with a National Instruments data acquisition 
board comprised the bulk of the data acquisition system 
used during the APU experiment. The acquisition board 
was a NB-MIO series board which featured an analog 
sampling rate of up to 100 KHz, 12-bit resolution, and 

16 digital Input/Output channels (ref. 46). LabVIEW®, an 
icon-based programming software tool was also installed 
onto the hard drive of the Macintosh computer to inter- 
face between the computer and the data acquisition 
board. Figure 15 is a How chart of the LabVIEW® data 
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acquisition and analysis routine as it was coded for the 
APU experiment. From the main routine, the experi- 
menter could take dynamic data, establish probe zero 
position, and adjust the traverse; other functions are listed 
in figure 15. An experimenter could move the traverse, 
monitor channels in very near real time by observing 
histograms and x-y plots, graphically observe the probe 
position in relation to the nozzle, and select nozzle type 
(circular or rectangular), all from the main menu of the 
LabVIEW® software program. 

2.3.12 On-the-Fly Plotting of Data- A personal 
computer (PC) was used to plot data on-the-fly during 
the APU experiment. Data were entered into a common 
spreadsheet program, Microsoft’s Excel® (ref. 47) as they 
were obtained. A back-up copy of the data was created by 
virtue of manually entering the measured temperature and 
pressure results. The primary data collection was com- 
pleted by the LabVIEW® data acquisition system. In 
addition to the spreadsheet itself, a current temperature 
plot and current pressure plot were each displayed on the 
PC monitor. Both the temperature and the pressure data 
plots were instantly updated upon data entry. By 
graphically visualizing the data as they were recorded, 
the experimenters gained additional insight into the 
validity of that data. For example, if the temperature and 
pressure decreased when it obviously should have 
increased, the experimenters would be alerted to check 
the probe’s position. 

2.3.13 Weather Station- Figure 16 is a photograph of 
the OARF’s weather station. The weather station, located 
approximately 150-ft (45.7-m) from the test site, provided 
continuous data such as ambient temperature, relative 
humidity, wind speed, and wind direction. Two three-cup 
anemometers and an aerovane can be seen mounted on 
the weather station tower in figure 16. The aerovane can 
measure both wind speed and wind direction; a three- 
blade propeller rotates at a rate proportional to the wind 
speed while the streamlined shape and vertical fin keep 
the blades facing into the wind (ref. 48). The digitally 
displayed results from the weather station were available 
in the OARF bunker. The ambient weather conditions 
were manually recorded at regular 30-min intervals. 

2.3.14 Agema Infrared Imaging System- An Agema 
Thermovision® Dual 880 Infrared Imaging System, 
shown in figure 17, was utilized to obtain infrared images 
of the APU exhaust plume. The imaging system employs 
two scanning infrared cameras that are similar in size to a 
hand-held video camcorder. The system is equipped with 
a data acquisition system, similar in size to a PC, and a 
color monitor and keyboard. Two sets of lenses are 
available for the infrared cameras; 20-deg lenses and 
7-dcg lenses provide a choice of field-of-view. The 


system operates on less than 400 W of 1 10- Vac power. 
The system is fully portable. The imaging system’s tripod 
is equipped with a fluid-filled rotating head to minimize 
jerky movement while scanning. 

The imaging system offers a variety of scanner modes; 
however, for the APU experiment the infrared cameras 
produced 25 images per second with 140 x 70 pixels per 
image. Each pixel of data is recorded with a full twelve- 
bit dynamic range. The infrared detector for the first 
camera is sensitive in the 2.5- to 5.5-jim wavelength 
band, which is referred to as the mid-waveband (MWB). 
The detector for the second infrared camera is sensitive 
from 4 to 13 Jim; however, it is predominately sensitive 
from 7 to 12 pm, which is referred to as the long- 
waveband (LWB). Figure 18(a) shows the spectral 
response function of the MWB infrared camera, and 
Figure 18(b) shows the response of the MWB camera 
when a narrow-band Filter has been added to reduce the 
spectral bandwidth. Figure 19 shows the spectral response 
of the LWB infrared camera. 

Each infrared camera employs a single cryogenically 
cooled detector element. The infrared detectors arc 
constructed of indium antimonidc (InAs) for the MWB 
camera, and mercury-cadmium-tclludde (MCT) for the 
LWB camera. InAs and MCT arc common infrared 
detector materials (ref. 49). The MWB detector provides 
a source temperature sensitivity of 0. 1 °C at 30.0 °C, 
while the LWB detector provides a sensitivity of 0.05 °C 
at 30.0 °C. Both scanners have an accuracy of ±2 percent. 
To keep these detectors cooled during continuous 
operation, liquid nitrogen (LN 2 ) must be added to the 
cameras on an hourly basis. 

The Agema system generates a multi-pixel infrared image 
by scanning a scene and digitally recording the infrared 
data to a hard disk (ref. 1 7). Up to 8 min of continuous 
MWB and LWB infrared data may be recorded on the 
750 Mbyte internal hard disk. A 1/4-in. tape drive 
provides a medium for secondary storage. Infrared 
radiation arriving from the scene is received at the 
detector element through a rotating-mirror assembly. An 
internal calibration is set on each scan line by sweeping 
over a LN 2 -cooled blackbody surface inside the camera. 
Measurement accuracy of the instrument is maintained by 
constantly monitoring internal temperature sensors, which 
compensate for detector-response drift by adjusting the 
gain. 

The color monitor provides a windows-based tool for 
displaying live or recorded images. Post-processing 
analysis tools such as spotmeters, isothermal lines, and 
image subtraction are available with the software required 
to run the system. Various functions provide a means of 
manipulating the images, i.e., adjustments to level and 
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range, so that the operator can enhance certain features of 
the How in the displayed image. 

2.3.15 Mylar Backdrop- A cold and uniform back- 
ground scene in an infrared image, such as those taken of 
the APU plume, can increase the quality of the image. For 
example, a cold or an apparently cold background will 
increase the contrast between the hot exhaust plume and 
the background. It is a goal that the low infrared emission 
of the background will not “wash-out,'” in the measured 
image, the small-scale structures of the exhaust plume, 
particularly those structures that exist near the shear layer 
created by the exhaust plume and the surrounding air. 

Cold and uniform backgrounds can be created in a variety 
of ways. A water-cooled backdrop is common. Typically 
anti-freeze is added to the water, and the water is 
refrigerated to maintain temperatures below freezing. 
However, one problem with water-cooled backdrops is 
condensation, which tends to form on the face of the 
backdrop when the relatively warm and moist ambient air 
comes into contact with the cold face. Condensation in 
the form of water droplets can destroy the uniformity of 
the backdrop. Also, water-cooled backdrops typically 
include a pump, a heat exchanger, and a water reservoir, 
all of which complicate transportation and portability. 

A reflective backdrop, shown in figure 20, was 
constructed for the APU experiment. In the visible 
spectrum as well as in the infrared, the reflective back- 
drop works like a mirror, reflecting the cold sky directly 
above the test site to the thermographic imaging system. 
The backdrop is relatively light-weight and it is portable; 
it can be easily moved by two people. The reflective 
material is aluminized Mylar, which is available at 
general hardware stores such as Tap Plastics in Mountain 
View, California (ref. 50). Mylar, a trademark of DuPont 
Corporation, is a thermoplastic film; aluminized Mylar 
has a thin film of aluminum on one surface of a Mylar 
sheet (ref. 51). The material is very thin, approximately 
0.003 in. (0.076 mm) thick. Tightly stretching the 
aluminized-Mylar over a wooden frame created a flat 
surface. Such a flat surface, tilted at an angle of 45 deg, 
is necessary to ensure that only the sky above the test site 
is viewed. Even small ripples in this mirror-like surface 
were to be avoided because they would reflect other 
scenes about the test site, including the hot nozzle. 
Reflecting only a limited portion of the sky was necessary 
to maintain a uniform background. 

The backdrop was positioned at a 45-deg angle with the 
ground plane for optimal reflection of the sky directly 
above the test site. In the MWB, the apparent temperature 
of the backdrop was approximately 55 °F (12.8 °C) and in 
the LWB the apparent temperature was -7 °F (-21.7 °C). 
Especially in the LWB case, the apparent temperature was 


lower, or better, than could be achieved by a comparable 
water-cooled system. Infrared images and further detailed 
results of the reflective-backdrop are given in Section 3. 

2.3.16 Midac System- A Midac spectrometer, utilizing a 
plane-mirror Michelson-type interferometer, was used 
one night during the axisymmetric APLI experiment to 
gather spectral data on the APU exhaust plume. The 
spectrometer is sensitive from 2.5 to 16 pm and has a 
resolution of 0.5 cm“’ . The system was used in con- 
junction with an infrared source and a laptop PC equipped 
with spectral analysis software (ref. 52). 

An informal agreement between NASA Ames Research 
Center and the Lawrence Livermore National Laboratory 
(LLNL) in Livermore, California, arranged for Charles 
Bennett of LLNL to visit NASA Ames for the purpose 
of obtaining spectral data on the APU exhaust plume 
(ref. 52). Bennett and experimenters at both NASA Ames 
and LLNL were optimistic that data collected from the 
spectrometer would be useful in identifying concentra- 
tions of exhaust plume constituents such as carbon 
dioxide and water vapor. With this objective in mind, 
Charles Bennett of LLNL visited Ames Research Center 
during the APU preliminary experiment in May 1994, to 
get a “feel” for the experiment. During the July 1994 lest, 
Bennett returned to Ames with the spectrometer to gather 
similar spectrometer data. Although results from the 
Midac spectrometer were inconclusive, a discussion of 
these results is given in Section 3. 

2.4 Experimental Procedures 

A typical night of experimentation would begin around 
1 1 :()() p.m. and end the following morning around 
8:30 a.m. The routine followed by the researchers during 
the 9.5-hr period was the same throughout each session. 
An outline of a typical test session is provided here to 
give the reader greater insight into how the test was 
conducted and how the data were obtained. For the 
axisymmetric plume survey, only four researchers were 
available to do the jobs described below. With only a 
few exceptions, which are explained below, the work load 
of completing these tasks was equally shared among the 
four workers. 

2.4.1 Open Gates to the OARF- The gates to the OARF 
were opened around 1 1:00 p.m. After the previous testing 
session, the gate to the facility was locked to discourage 
passers-by from interfering with equipment which 
remained set up for the duration of the test. 

2.4.2 Open Trailer and OARF Bunker- Next, the 
trailer and OARF bunker were unlocked and the spot- 
lights were turned on. The OARF has a large underground 
support bunker which is used during larger tests; the 
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bunker was constructed underground so that no structural 
interference would occur between the building and any 
above-ground test. This facility is heated and is equipped 
with a microwave oven and a rest room. Also, the digital 
readout of data from the weather station is available in 
this bunker. Two large mercury spotlights are powered 
from this room; the spotlights served to adequately light 
the test rig. 

2.4*3 Uncover Test Equipment- After the mercury 
spotlights were illuminated, the canvas and plastic 
tarpaulins were removed from the traverse rigs and other 
equipment located outdoors. The tarpaulins were used to 
cover equipment such as the traverse rigs, the power 
rectifier, and the motors and motor controllers to protect 
these items from rain and dust. 

2.4.4 Set Up the Aluminized Mylar Backdrop- The 

next task required two people to attach the aluminized 
Mylar sheet to the backdrop structure. The aluminized 
Mylar sheet is very thin and consequently cannot 
withstand the heavy breezes which blow during the 
afternoons between test periods. For this reason, after 
each morning of testing the Mylar was removed from the 
support structure. In order for the Mylar to reflect the cold 
sky as well as possible, it was necessary to tightly stretch 
the Mylar over the backdrop support to remove any 
wrinkles. The Mylar was held in place by duct tape which 
was spaced at regular intervals along the perimeter of the 
backdrop; one person would stretch and a second person 
would secure the Mylar to the backdrop with tape. 
Because the infrared emission from the duct tape was 
much greater than the radiation reflected from the 
backdrop, the duct tape, placed on top of the Mylar, 
provided a sharp contrast in the infrared. Consequently, 
the regular spacing of the duct tape provided a method for 
scaling the plume in the infrared image; the known and 
regular spacing of the duct tape resembles the markings 
on a ruler against the Mylar backdrop in the infrared 
images. This process of attaching the Mylar to the 
backdrop structure required at least twenty minutes. 

2.4.5 Attach Probe to the Three-Dimensional 
Traverse- The next step involved attaching the 
combination temperature and pressure probe to the 
mounting block on the small three-dimensional traverse. 
The combination probe is a delicate and expensive 
transducer and was consequently treated with great care. 
The probe is mounted to a large aluminum block, which 
in turn is mounted to the small traverse. Two tightening 
screws secure the rear one-third of the probe to the 
aluminum block. The large aluminum block adds rigidity 
to the probe, helping the probe hold its position within the 
quasi-steady flow. Once the probe had been securely 
fastened to the aluminum block, the typc-K thermocouple 


wire and pressure tube leading from the probe were 
connected to the wire and tube leading to the SCU. At 
this time a second probe was attached to the large traverse 
rig to measure ambient conditions. 

2.4.6 Reattach Motor-Controllers- Next, the motor 
controller power cables were reattached and the con- 
trollers were turned on. The power cords were removed 
daily from the motor controllers in an effort to avoid 
accidents. A digital display on the face of the controller 
confirmed successful power-up self-test diagnostics. 

2.4.7 Set Up the Infrared Imaging System- Each night 
the same person was responsible for setting up the 
infrared imaging system. These responsibilities included 
assembling the camera system, attaching the appropriate 
cables and lenses, obtaining the liquid nitrogen and 
attaching the cameras to the sturdy tripod. It was impor- 
tant that the subtleties of camera location and height did 
not vary from night to night. For this reason, duct tape 
was pasted on the concrete pad to mark the exact position 
of the legs of the tripod. A level on the head of the tripod 
was also used to establish camera position. 

It was anticipated that by using the same operator to do 
these jobs, the job would be done consistently from night 
to night. This is especially true when taking data. When to 
take data, the duration of a data run, and the selection of 
lenses are each important operator decisions. By main- 
taining a consistent procedure throughout the test, the best 
possible data set was obtained. Also, post processing was 
made easier by consistent data filename selection. Once 
the equipment had been set up, liquid nitrogen was added 
to the dewars in each of the cameras. After allowing a few 
minutes for the infrared detectors to cool, the system was 
turned on. By comparing the infrared scene created by the 
cameras to a scene taken from a previous night, final 
position adjustments were made to the camera to ensure 
consistent results from one testing session to the next. 
Finally, the experimenter would check the available space 
on the hard disk of the Agema system to make sure that 
sufficient data storage space was available for the 
evening’s scheduled tests. 

2.4.8 Power-up Computer Systems- An experimenter 
would turn on the power to the PC and Macintosh 
systems in the trailer. Successful boot-up of Microsoft 
Excel® and the National Instruments LabVIEW® 
software was verified. 

2.4.9 Take Weather Data- At this point, the 
experimenters would begin taking weather data on the 
half-hour. The digitally displayed weather data such as 
wind-speed, wind-direction, ambient temperature, and 
relative humidity were available in the OARF bunker. An 
experimenter would walk into the bunker every half-hour 
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and record the displayed values on a pad of paper. The 
first data point was usually taken around 1 1 :30 p.m., and 
the last point taken around 7:30 a.m. 

2.4.10 Pre-Experiment Test Meeting- A pre-experiment 
test meeting was held each night to discuss the schedule 
of events. Personnel safety information as well as 
experiment objectives were reviewed. Results Irom the 
previous day’s testing were shown and posted on a cork 
board inside of the trailer. 

2.4.11 Check-Out of Data Acquisition System- 

Because of its important function, a thorough check of 
the Macintosh computer system and LabVIEW® data 
acquisition system was performed each night before 
testing began. One key item on the checklist was the 
verification that adequate data storage space was available 
on the removable hard disks that served as the primary 
data storage medium. On one particular night during the 
axisymmetric test, several hours of data were lost because 
we were writing to a full disk and had not programmed 
LabVIEW® to warn us of this error. Luckily, the average 
temperature and pressure data points were saved because 
of the manual data entry into the Excel® spreadsheet. 
Another important detail not to be overlooked was probe 
positioning. Occasionally, the probe would erroneously 
not move when it should have moved. For this reason, a 
complete checkout of the LabVIEW® system involved a 
visual inspection of the probe’s position. 

2.4.12 Turn On Fuel Pumps- A long fuel hose, 
approximately 100-ft (30.5-m) long, delivered fuel from a 
large, 10, 000-gal (37,854 I) underground fuel tank to the 
APU. An experimenter would turn on the fuel pump at 
the pumping station and leave the fuel valve and pump 
switch accessible as a precautionary measure. Visual 
inspection of an in-line fuel pressure gauge was required 
to make certain that the line gage pressure of 12 Ibf /in 2 . 
(82.7 kPa) remained constant throughout the experiment. 

2.4.13 Turn On 28-Vdc Power Rectifier- A 480- Vac - 
to-28-Vdc power rectifier was used to convert the OARF 
480-Vac electric service to the 28 Vdc power required to 
operate the APU. The alternative power source, a stand- 
alone 28-Vdc power supply, was extremely noisy 
(=120 dB); the rectifier provided a relatively quiet 
alternative to this option. A 28-Vdc power source was 
required to run the fuel pump on the APU. For this 
reason, the rectifier was turned on before the APU was 
started. 

2.4.14 Warm-up of APU- Next, an operator would start 
the APU and allow it to warm up. A visual inspection of 
the revolutions-per-minute (RPM) gauge and the exhaust- 
gas temperature (EOT) gauge was required to make 
certain the APU started correctly. A kill-engine switch. 


located near the RPM and EGT gauges, was actuated 
in the event the APU did not start correctly. Bad starts 
were sometimes caused by pockets of air in the fuel line. 
After the APU had “warmed-up,” approximately fifteen 
minutes after start, the APU was turned off. At this point 
a relatively stable operating condition had been estab- 
lished. Also, the stainless steel exhaust nozzle and flow 
straightening assembly had heated to an equilibrium 
temperature and consequently expanded. 

2.4.15 Level and Straighten Exhaust Nozzle- The level 
of the traverse was rechecked and if necessary, the 
exhaust nozzle of the APU was leveled and straightened. 
The level of the traverse rig was checked; it usually did 
not require adjustment after the initial leveling. Also, the 
level of the exhaust nozzle was checked. The straightness 
of the nozzle was checked using a plumb, which was 
hung from the probe and traverse rig. Making sure the 
nozzle was straight and level, as well as assuring that the 
probe traversed level and straight, was an important issue 
in successfully surveying the plume. 

2.4.16 Establish Probe “Home” Position- Alter the 
level of the traverse rig and exhaust nozzle had been set, 
the probe was positioned at a “home” position. The 
position was measured near the exhaust nozzle exit. If the 
probe was positioned near this exit before the nozzle had 
expanded, the probe could be damaged by expansion of 
the nozzle. The “home” position was verified by visual 
inspection using a steel ruler. 

2.4.17 Re-Start APU and Allow for “Warm-Up”- 

Next, the APU was restarted and thermal equilibrium 
attained. By viewing the on-line probe temperature 
reading, at this point centered at the exhaust exit, the 
experimenter could tell when a steady temperature had 
been obtained. Once the steady temperature was reached, 
five minutes elapsed before any data were recorded. The 
five minutes allowed for further nozzle expansion. At this 
point the APU was left running until some event 
warranted its shut down. 

2.4.18 Take Data Point “Home” Position- The first data 
point recorded during a testing period usually occurred at 
the “home” position. This event served as a final check o( 
successful system operation and gave the experimenters 
an indication of how the APU was running for that 
session. Depending upon the ambient temperature, the 
temperature at the “home” condition could vary as much 
as 20 °F (1 1 K) from session to session. 

2.4.19 Take Infrared Data of the APU Exhaust Plume- 

Alt hough one entire shift was dedicated to taking intrared 
data of the APU exhaust plume, additional infrared data 
were taken throughout the APU test. These data included, 
for example, images of the backdrop, start-up and 



shut-down images, and images taken with and without the 
probe in the flow. 

2.4.20 Begin Taking Data- Beginning with the first 
measured data point, the experiment was underway. The 
interaction between the Lab VIEW® operator and the 
Excel® operator played an important role in the data 
acquisition process. To begin, the data points to be 
obtained during the nights work were listed on the Excel® 
spreadsheet. One column of that spreadsheet listed the 
axis movements required to reach the next data point. 

For example, that column may have read move +0.25 in. 
(+0.635 cm) in the x-direction. The Excel® operator 
would announce this traverse movement to the 
LabVIEW® operator, who in turn would click the 
appropriate icon on the terminal of the Macintosh to 
move the probe to that position. After a few seconds, an 
audible tone would alert both operators to the fact that the 
probe was on position. The current x-y-z position was 
displayed on the LabVIEW® screen so that the Excel® 
operator could double-check the probe position with the 
position listed on the spreadsheet. 

Next, the LabVIEW® operator would watch the displayed 
average temperature. When that temperature had 
stabilized, i.e., when it was no longer increasing if going 
into the plume, or vice-versa for exiting the plume, then it 
was time to take data. The LabVIEW® operator would 
click the “take data” icon at this time. One hundred 
consecutive temperature and pressure measurements were 
then written to a removable hard disk. The average 
temperature and pressure measurements were displayed 
on the Macintosh terminal, where the average value was 
computed using the one hundred measured values. 
Extraneous data were removed by clipping those data 
that fell outside a two-standard-deviation window. The 
LabVIEW® operator would then announce those average 
values to the Excel® operator, who would in turn enter 
those values into the spreadsheet. Next, the Excel® 
operator would announce the subsequent traverse move- 
ment. This process was repeated until the experiment was 
stopped for some reason. 

2.4.21 Visual Inspection of the Probe Position- 

Becausc no feedback was available on the position of the 
moving probe, a periodic visual inspection of the probe 
position was required. Two windows in the trailer gave 
the experimenters seated at the computers a view of the 
probe. However, it was generally required to look at the 
probe from close range to verify that it was correctly 
positioned. 

2.4.22 Save Data Regularly- Because the standard 
timed-backup routine for the Excel® spreadsheet saved 
data too infrequently to satisfy the test crew, the experi- 
menters frequently reminded each other to save data. 


Enough data was collected and entered into the spread- 
sheet between the standard timed-backups to assure that 
experimenters eagerly took a more active role in saving 
the data. Data were saved to a variety of filenames as an 
added precaution. 

2.4.23 Turn Off the APU- When the final data point had 
been recorded, the APU was shut down. During a testing 
session, the APU was sometimes shut down for various 
reasons that included, but were not limited to, problems 
such as incorrect probe positioning, insufficient data 
storage space, high winds, or problems with the data 
acquisition system. Otherwise, the APU was left on until 
the end of the testing session. 

2.4.24 Record Position of the Test Equipment- In an 

effort to ensure the repeatability of the experiment from 
one testing session to the next, daily measurements of the 
equipment position were made. For example, the position 
of the APU, the backdrop, the camera system and the 
traverse were all verified by using a tape measure. 
Photographs taken during the experiment also helped to 
document equipment position. 

2.4.25 Turn Off the 28-Vdc Power Rectifier- After the 
APU had been shut down, the 28-Vdc power rectifier was 
turned off by actuating the appropriate switch on the 
instrument panel. The relatively hot rectifier had to cool 
down before it could be covered with a tarpaulin. The 
power cable leading from the rectifier to the APU was 
rolled up and safely stored on a rack on the side of the 
rectifier. The power cable leading from the rectifier to the 
480- Vac power source was likewise stored in a rack on 
the side of the rectifier. 

2.4.26 Secure the Fuel Pumping Station- After each 
testing session, the fuel pump power was turned off at the 
fuel pumping station. Next, the cover was returned to the 
fuel pit area which leads to the flow valves and fuel 
pump. This cover was removed during pump operation to 
give safety personnel easy access to the fuel pit area in the 
event of an accident. Finally, the fuel pump power switch 
was secured by locking the face-plate leading to the 
switch. 

2.4.27 Dismantle the Aluminized Mylar Backdrop- 

The thin aluminized Mylar sheet was removed from the 
backdrop and stored in the trailer after each testing 
session. Care was taken to roll up the Mylar sheet only 
when the sheet was completely dry. Moisture, caused by 
an infrequent morning dew, would cause the Mylar to 
stick to itself when the sheet was rolled up, creating the 
possibility that the aluminum would separate from the 
Mylar when the sheet was unrolled. The duct tape used to 
secure the Mylar to the frame was carefully removed to 
prevent the sheet from tearing. The aluminum removed 
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by the duct tape posed no threat to the function of the 
backdrop because the duck tape was fastened to the Mylar 
only on the edges of the Mylar sheet. The large wooden 
backdrop frame was set on its side to prevent accidental 
damage to the frame caused by the strong winds that blew 
during the afternoons between tests. 

2.4.28 Store the Infrared Imaging System- The 

infrared imaging system is an extremely delicate and 
expensive instrument. One cost estimate for replacing the 
system exceeds 200,000 U.S. dollars (ref. 54). For this 
reason, after each testing session the system was carefully 
repackaged and stored in a different building. 

2.4.29 Make Back-Up Copies of Survey Data- The 

complete temperature and pressure data sets were saved 
on several removable optical diskettes. The average 
temperature and pressure data set was saved onto the hard 
drive of the PC. Back-up copies of this reduced data were 
made after each testing session by copying those files to a 
diskette. The diskettes were stored at a different location 
between tests. After the back-up copies were made, the 
computers were turned off. The keyboards to the 
computers were removed and stored in a different 
building for added security. 

2.4.30 Secure OARF Bunker- The underground control 
room at the OARF was tidied and secured. The mercury 
spotlights were turned off. The last weather data were 
recorded and the weather conditions for the session were 
cataloged. 

2.4.31 Secure the Facility for Subsequent Tests- The 

final wrap-up of events began by securing loose items in 
the trailer. The test equipment on the concrete pad was 
covered by replacing the plastic and canvas tarpaulins. 
Duct tape and nylon ropes were used to secure the 
tarpaulins. Around 8:30 a.m., the personnel would depart 
the OARF. The gate to the OARF was locked at this time. 

3. Results 

This section presents the results from the axisymmetric 
APU exhaust plume survey. Temperature and pressure 
plots, weather data, infrared data, and spectrometer data 
arc shown and discussed. 

3.1 Choice of the Coordinate System 

A right-handed Cartesian coordinate system was selected 
for describing the APU exhaust plume. The x,y-axis plane 
is parallel with the exit plane of the exhaust nozzle. The 
positive z-axis begins at the exit plane of the nozzle and 
proceeds downstream along the centerline of the exhaust 


flow. Figure 21 shows the coordinate system in relation to 
the axisymmetric nozzle. 

3.2 Summary of Temperature and Pressure Data 
Taken 

Recall that two experiments were conducted (see 
Section 2.1 ) in which the How exiting the axisymmetric 
nozzle was surveyed. The first experiment was a pre- 
liminary experiment conducted in May 1993, and the 
second experiment was a dedicated axisymmetric nozzle 
survey conducted in July 1993. Table 1 summarizes the 
data obtained during the preliminary experiment and 
tables 2 and 3 summarize the data obtained during the 
dedicated test. Notice from these tables that several 
repeatability runs were conducted to establish confidence 
in the measured data. 

3.3 Uncertainty of Temperature Measurements 

The temperature of the APU exhaust plume was measured 
using the combination temperature and pressure probe 
described in Section 2.3.7. The temperature probe was a 
type-K thermocouple, made of Chromcl-Alumel. The 
probe was calibrated using the Omega Hand-Held 
Calibration Unit, Model CL23-JKT Calibrator (rel. 44), 
which is described in Section 2.3.8. The calibration unit 
provided an accuracy of 0. 1 percent of the recorded 
temperature plus 1 °F. The thermocouple was calibrated 
one week prior to both the preliminary experiment in May 
and the dedicated APU experiment in July. 

3.4 Uncertainty of Pressure Measurements 

The pressure transducer portion of the combination 
temperature and pressure probe is also described in 
Section 2.3.7. The Kiel probe calibration unit is discussed 
in Section 2.3.9. The calibration unit provided an 
accuracy of 0.05 percent of full scale (30 psid max.), 
which implies accuracies of 0.015 psid. Table 4 lists a 
comparison between the pressures entered into the 
calibration unit and the corresponding pressures indicated 
by the Kiel probe. The percent difference between the 
pressure indicated by the Kiel probe and the pressure 
maintained by the calibration unit is listed in the third 
column of the table. The majority of the values listed in 
the third column show agreement within ±0.5 percent. In 
addition, it should be reported that maintaining a stable 
pressure on the calibration unit proved to be a fairly 
difficult task; the pressure would fluctuate between 
±0.002 psid. This value is, however, within the implied 
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accuracy of 0.0 1 5 psid. The calibration unit itself was 
calibrated just prior to these measurements. 

3.5 Data Clipping and Data Averaging 

As mentioned in Section 2.3.1 1, the Labview® data 
acquisition system played a key role in the experiment. 
When the experimenters obtained a data point, one 
hundred essentially-instantaneous temperature and 
pressure measurements were recorded. These data were 
reduced to a single average temperature and average 
pressure value. The average temperature and pressure 
values were computed by averaging those data that fell 
within a two-standard-deviation, 2-a, window. The 
standard deviation, a, of the numbers was computed 
using the formula 



In equation 1 , n is the number of data points and x is 
either the temperature or pressure. Once a was computed, 
those data that fell outside of a 2-a window were dis- 
carded. The process of removing extraneous data in this 
fashion is standard practice among experimental scientists 
at NASA Ames (ref. 36). The average temperature or 
pressure value, x , was computed by averaging the 
remaining data by 



where n* is the number of remaining data points. 

An illustration of the clipping process is shown in 
figure 22. Here, data falling above 745.2 °F or below 
742.9 °F have been “clipped” or removed. The variance 
of these data, temperature or pressure, is a consequence of 
taking measurements on a dynamic system. Although the 
flow was straightened and the system was relatively 
stable, it is reasonable to expect some fluctuation in the 
temperature and pressure values. This phenomenon was 
also observed in the infrared data and is explained further 
in Section 3. 10. 1 . 

3.6 Probe Heating Study 

A simple experiment was conducted to investigate the 
effect of probe healing, i.e., to determine the effect, if 
any, the heating of the probe by the exhaust stream had on 
the measured temperature and pressure data. For this test, 
the probe was caused to traverse the flow by inserting the 
probe into the How from the ambient air to the centerline 


of the exhaust stream. Temperature and pressure 
measurements were obtained at regular intervals during 
the traverse. After reaching the center of the flow, the 
probe was traversed backwards along the insertion line, 
i.e., the probe was traversed from the center of the flow to 
the ambient air. Temperature and pressure measurements 
were made at the same locations along the traverse line. A 
comparison of the data between the two traverses should 
reveal the effect heating of the probe by the exhaust 
stream would have on the measured data. For example, 
for the same spatial position in the exhaust stream, the 
probe would likely yield a cooler temperature on the 
entering traverse than on the exiting traverse due to 
heating of the probe by the exhaust stream. 

Figure 23 shows the results from this experiment obtained 
during the rectangular nozzle case (no such study was 
conducted for the axisymmetric case). The rectangular 
nozzle is different from the axisymmetric nozzle in that 
less of the probe length is heated by the flow. However, 
an approximately ten percent warmer centerline tempera- 
ture present in the rectangular nozzle should exaggerate 
the effects of heating. Figure 23(a) shows the measured 
temperature data during the entering and exiting traverses, 
while figure 23(b) shows the measured pressure data. The 
data obtained during the entering and exiting sweeps 
agree remarkably well. The favorable comparison 
between these data gives us confidence that the data 
sampling rate was sufficiently slow enough to allow the 
probe to properly cool on the exiting traverse. 

3.7 Temperature and Pressure Data 

The results from the temperature and pressure 
measurements arc shown in figures 24-35. Figure 24 
shows results at the exit plane, figure 25 shows results at 
one-half diameter, and figures 26-35 show results from 
one to ten diameters downstream at regular one diameter 
intervals. Where applicable, data from the preliminary 
experiment are shown. In general, the y- or x-axis 
location is shown on the horizontal axis of each graph, 
and the temperature or pressure is shown on the vertical 
axis. Each figure is a summary of all average temperature 
and pressure measurements obtained at a particular z-axis 
location. Recall that a summary of these data is given in 
tables 1-3. Appendix A includes forty tables that list the 
experimentally obtained temperature and pressure data 
that are presented in figure 24 -35. These tables are 
referenced in the caption of each figure. 

3.7.1 Preliminary Experiment Data Compared to 
Dedicated Experiment Data- In figure 24, temperature 
and pressure data are presented for both y-axis and x-axis 
sweeps. These data are compared to measurements 
obtained during the preliminary experiment and are 
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consequently labeled, for example, “prelim, run 1.” 

Notice that the trends of the curves from both the pre- 
liminary experiment and the dedicated experiment are 
similar. Because Run 1 and Run 2 fall nearly on top of 
one another, we have an additional degree of confidence 
in the measured data. The same is true for the preliminary 
experiment data. 

3.7.2 Higher Temperatures and Pressures Observed 
in the Dedicated Experiment- A comparison between 
the preliminary experimental data and data from the 
dedicated test reveals that higher temperatures and 
pressures were obtained during the dedicated test. This 
phenomenon can be observed in figures 24-30. It may be 
hypothesized that the APU ran warmer in the dedicated 
experiment due to longer run times and a higher fuel line 
pressure. During the preliminary experiment fuel was 
delivered to the APU from a small 10 gal (37.8 1) fuel 
tank which is shown in figure 8. Because of the small 
capacity of this fuel tank, the APU had to be stopped 
approximately every forty minutes for refueling, allowing 
the APU to “cool off." Recall that during the dedicated 
APU experiment, fuel was delivered to the APU from a 
large fuel tank, 10,000 gal (37854 1), and fuel pumping 
station (see Section 2.4.12). Because the APU did not 
have to be shut off for refueling reasons during the 
dedicated experiment, the APU was left running for 
longer periods of time. It is hypothesized that the warmer 
exhaust temperatures and higher exhaust pressures that 
are observed in the dedicated experiment are a conse- 
quence of the longer run times. 

3.7.3 Asymmetries Observed in the Axisymmetric 

Plume- In tracing the measured pressure values from 
the ambient air to the centerline of the exhaust plume, 
a monotone increasing line would be indicative ot a 
symmetric flow. Likewise, data obtained on the x-axis 
should compare with data obtained on the y-axis if 
the APU exhaust plume were truly symmetric. Unfortu- 
nately, a hump in the pressure curve is observed in 
figures 24(a)-28(a). Likewise, an asymmetry is observed 
in figures 24(c)-28(c). This asymmetry is likely due to 
some error in the construction of the tube bank, a buildup 
of soot on a straightening screen, or a disturbance induced 
in the flow downstream ot a welded seam. These asym- 
metries tend eventually to wash-out by four diameters 
downstream, as may be observed in figures 29-35. 

3.7.4 Limitation of x-Axis Data to Six Diameters 
Downstream- Both x- and y-axis data arc shown in 
figures 24-31, which show those graphs corresponding to 
the exit-plane data through data obtained at six diameters 
downstream of the exit. Beyond six diameters, that is 
from seven diameters downstream to ten diameters 
downstream, only y-axis traverses are presented. For the 


corresponding figures 32-35, the x-axis data were not 
measured due to a limitation on the freedom of movement 
of the small traverse in that direction. 


3.8 Velocity Profile 

The velocity at point «, u n (ft/s), is given by 


u n = K P Cp (^p) 
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In equation 3, Kp is the Pitot tube constant given by 
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Cp is the Pitot tube coefficient (dimensionless), Ap is the 
velocity head measured by the probe (in. H 2 O), T n is the 
temperature (°R), P n is the pressure (mm Hg), which is 
computed by 

P n = P a t m + A p ( 5 ) 

where P atm is the atmospheric pressure, and M is the 
molecular weight of the gas (Ib/lbmol). 

By using equation 3 to compute the velocity, the 
temperature and pressure profiles presented in Section 3.7 
can be converted to velocity profiles. Figures 36-41 show 
the velocity profiles from the exit plane to ten diameters 
downstream. The y-axis profiles were rotated about the 
centerline of the exhaust plume because, unlike the x-axis 
traverses, the y-axis traverses were taken at each survey 
plane from the nozzle exit to ten diameters downstream. 

Figure 36(a) shows the velocity profile at the exit plane. 
Similarly figure 36(b) shows the velocity profile at one- 
half diameter downstream. Both figure 36(a) and 
figure 36(b) show steep profiles. The dip in the velocity 
profile at the centerline is a consequence of the pressure 
nonuniformity discussed in Section 3.7.3. 

Figure 37(a) and figure 37(b) show velocity profiles 
for one and two diameters downstream, respectively. 
Figure 38(a) and figure 38(b) show similar profiles for 
three and four diameters downstream, respectively. 
Between figure 37(a) and figure 38(b), evidence that the 
exhaust plume is spreading as it progresses downstream is 
clearly visible. 

As described in Section 3.7.3, the effect of the asymmetry 
begins to wash out beginning at four or five diameters 
downstream. Evidence of this occurrence is apparent in 
figures 38(b), 39(a), and 39(b), the figures corresponding 
to the velocity profiles at four, five, and six diameters 
downstream, respectively. 
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Figures 40 and 4 1 display the velocity profile at nine and 
ten diameters downstream, respectively. The peak 
velocity at ten diameters downstream is 250 ft/s 
(76.2 m/s). This velocity is slightly less than one-half of 
the peak velocity at the exit plane of 530 ft/s (161.5 m/s). 
This peak velocity, shown in figure 36(a), corresponds to 
a Mach number of 0.47. 

3.9 Weather Data 

As described in Section 2.3. 1 3, the weather station at the 
OARF was used to gather data such as ambient tempera- 
ture, atmospheric pressure, relative humidity, wind speed 
and wind direction. Tables B1 through B5, given in 
Appendix B, list these data for July 12, 1993, to July 16, 
1993, respectively. The ambient weather conditions are 
listed for approximately every half-hour of the dedicated 
axisymmetric experiment. The time-averaged values of 
those parameters are computed and listed for each day. 

The time-averaged ambient temperature over the five-day 
experiment varied between 57 and 58 °F (13.8 and 
14.4 °C). In general, the temperature dropped through a 
testing session, reaching its lowest value at daybreak. The 
relative humidity varied between 80 and 95 percent, in 
general increasing through the night. The wind speed 
varied between 4 and 6 mi/hr (6.4 to 9.7 km/hr). The 
consistent ambient temperature, a stable relative 
humidity, and a low wind speed are desirable conditions 
for this type of outdoor experiment. Operating at night 
between 12 midnight and 8 a.m. allowed us to take 
advantage of these conditions. 

3.10 Infrared Image Data 

A large selection of infrared image data was recorded 
during the axisymmetric APU exhaust plume survey 
experiment. Appendix C includes a list of these data Hies. 
Each table in Appendix C lists the infrared data that were 
measured during a given night of the dedicated APU test. 
Over 1 30 data files were recorded over the duration of the 
experiment. The images include both MWB and LWB 
data. Data files include background images, images from 
different view points, and images measured with and 
without the nozzle shield. These data files were recorded 
using the 20-deg lenses. The infrared images that are 
presented and discussed in this section of the report 
represent only a small portion of the actual data that has 
been obtained. However, the data shown are representa- 
tive of what is stored at NASA Ames Research Center on 
magnetic media. Those readers who have a need for an 
electronic version of these data may obtain them as well 
as the temperature and pressure data, by contacting the 


appropriate civil servant at NASA Ames, as discussed in 
Section 5. 

3.10.1 Field Versus Frame- The Agema 880 Infrared 
Imaging System creates an infrared image by scanning a 
scene with a rotating-mirror and a fixed-mirror assembly 
embedded within the camera. The system is capable of a 
variety of scanning rates; however, only one scanner rate 
was employed for this experiment. Specifically, the 
scanners produced 6.25 images per second with 

140 x 280 pixels per image. One image is referred to as 
one “frame” of data, where each “frame” is comprised of 
four interlaced “fields” of data. As shown in figure 42, 
one frame of data ( 140 x 280 pixels) is constructed by 
interlacing together four fields of data (140 x 70 pixels). 
The interlacing process creates the first four seventy-pixel 
rows of the frame, by using the first rows of the first, 
second, third, and fourth measured fields of data. The 
next four seventy-pixel rows of the frame arc created by 
using the second row of the first, second, third, and fourth 
measured fields. This process is continued until the 
140- x 280-pixel image is created. 

3.10.2 Measurement of a Dynamic System- Because 
each field of data takes 0.04 s to obtain, the relatively 
short time span allows one the option of capturing the 
dynamic nature of a flow field such as the APU. 

Figure 43 shows three consecutively measured fields of 
filtered MWB infrared data. Each image of the APU 
exhaust plume shown in figure 43 looks slightly different. 
The difference between these images illustrates an 
important concept: while many precautions were taken 

to ensure that steady How conditions were achieved, the 
APU exhaust plume is still a dynamic, or time-varying 
system. 

3.10.3 MWB Infrared Data- The three infrared images 
shown in figure 44 are MWB data. Figure 44(a) is a 
MWB field of data. Figure 44(b) is a MWB frame of data 
that has been created by interlacing four images like the 
one shown in figure 44(a). The MWB image shown in 
figure 44(c) is the average of six frames of data. The 
horizontal zigzag lines, or cuts, visible on the outer 
fringes of the exhaust plume in figure 44(b) and 

figure 44(c) are a consequence of building the image by 
interlacing together several frames of data. An artifact of 
imaging adynamic system with a scanning instrument, 
these horizontal lines tend to disappear as more fields of 
data are averaged together; see Section 3.10.6. Comparing 
figure 44(b) with figure 44(c), these horizontal lines arc 
seen to nearly disappear. 

3.10.4 Filtered MWB Image Data- The infrared images 
in figure 45 are filtered MWB data. Figure 45(a) shows 
one field of data, figure 45(b) shows one frame of data, 
and figure 45(c) shows an image averaged over six frames 
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of data. The color scale shown at the bottom of the figure 
applies to all three images. Recall that the filtered 
response of the MWB camera is shown in figure 18(b). 

One goal of this report is to provide a data base for 
validating Turk’s (ref. 25) infrared image prediction 
algorithm described in Section 1.2. The filtered MWB 
image shown in figure 45(c) represents an image ideal for 
this purpose. In figure 45(c) the backdrop is uniform, the 
spread angle is visible (see Section 3.10.8), and the probe 
has been positioned outside of the exhaust stream. An 
intensity scale has been added to this image. This image 
is compared in Section 4 to an infrared image predicted 
using Turk’s algorithm. 

3.10.5 LWB Image Data- Figure 46 shows LWB images 
of the AFU exhaust plume. Figure 46(a) shows one field 
of data, figure 46(b) shows one frame of data, and 
figure 46(c) shows an infrared image averaged over six 
frames of data. The wavy streak visible at the top right- 
hand corner of each of these three images is the edge of 
the aluminized-Mylar backdrop. The edge is visible in 
these LWB images because the cameras were positioned 
to get the best possible plume and backdrop alignment for 
the MWB camera. Consequently, the LWB image shows 
some of the backdrop. Despite the appearance of the 
backdrop on the edge of the LWB images, the backdrop is 
located directly behind the AFU plume. The color scale 
shown at the bottom of figure 46 applies to each of these 
three images. 

3.10.6 Advantages of Time-Averaging- Figure 47 
shows two MWB infrared images, each averaged using 
five hundred frames of data. This represents an average 
of data over eighty seconds in time. The horizontal cuts 
described in the preceding sections are less apparent in 
these images, making the plume appear smooth with less 
dynamic fluctuations. The combination temperature and 
pressure probe is visible as a bright point source in 
figure 47(a). The probe is obviously being heated by the 
flow. Visible above the probe is another “hot spot.” This 
second “hot spot” is a reflection coming from the under- 
side of the unrolled aluminized Mylar tube. 

An additional point of interest about figure 47(a) concerns 
the white, or hot, area near the nozzle exit. For most of 
the infrared images obtained during the APU experiment, 
a foreground shield, or reflective barrier, was placed 
between the APU nozzle and the camera. Figure 47(b) 
shows a MWB infrared image where the foreground 
shield has been used to obscure the APU nozzle. Hot 
parts such as metal nozzles tend to radiate with a much 
larger intensity than do exhaust plumes. For this reason a 
shield is placed between the nozzle and the camera to 
prevent the additional energy from the nozzle from 
reaching the detector of the camera and potentially 


saturating the image. This has clearly occurred in 
figure 47(a). With the exception of figure 47(a), which is 
presented for illustration only, each of the infrared images 
presented in this section were taken with the foreground 
shield in place. 

3.10.7 Infrared Image of the Backdrop- The alumi- 
nized Mylar backdrop, as discussed in Section 2.3. 15, 
provides a uniform, relatively cold, and substantially 
large area for imaging of the APU exhaust plume. 

Figure 48(a) shows a MWB infrared image of the 
backdrop where the nozzle foreground shield is visible 
at the lower portion of the image. Figure 48(b) show's a 
LWB infrared image of the aluminized Mylar backdrop. 
Visible on either side of the backdrop is the horizon. 

The difference between the relatively uniform radiance 
(reflected) of the backdrop and the relatively non-uniform 
radiance of the horizon dramatically illustrates the 
advantages of using a backdrop. The average apparent 
temperature of the backdrop shown in figure 48(a) is 

57 °F (13.9 °C) with a standard deviation of 5 °F. The 
corresponding apparent temperature of the backdrop in 
the LWB is -7 °F (-21.7 °C) with a standard deviation of 
2 °F. Each of these calculations were computed using 
approximately twenty-six thousand pixels. The apparent 
temperature is lower in the LWB than in the MWB 
because the LWB radiation comes from higher in the 
atmosphere due to less absorption in that band by 
atmospheric gases. 

Because the aluminized-Mylar backdrop works by 
reflection, clouds and airplanes that pass above the test 
scene can affect uniformity and apparent temperature. 

For example, an overhead cloud can cause the apparent 
temperature to fluctuate as much as 30 °F. For this reason 
the backdrop worked best on cold and cloudless nights. 

3.10.8 Plume Spread Angle- The infrared images shown 
in figures 44-48 illustrate the plume spreading angle. 
Hossain and Rodi (ref. 55), Chen and Rodi (ref. 56), and 
Bejan (ref. 57) predict spread angles between 7.5 deg and 
9 deg for turbulent, buoyant, vertical jets. The spread 
angle of this jet is approximately 7 deg from vertical. This 
result is consistent with Hardman’s observations (ref. 20). 

3.11 Validity of the Measured Infrared Image Data 

The infrared images presented in this report are accurate 
only to within the limits of the Agema 880 Dual Band 
scanning infrared camera system. The MWB camera is 
quoted by Agema as having a source sensitivity of 0. 1 °C 
at 30.0 °C. The LWB camera has a source sensitivity of 
0.05 °C at 30.0 °C. The scanners have an accuracy of 
±2 percent. Keep in mind that these data are quoted for 
a new camera system and that detector elements can 
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degrade with use. The Agema system used in the APU 
experiment was recalibrated by the Agema company just 
prior to the dedicated APU exhaust plume survey. 

Proper use and storage of delicate instruments such as the 
Agema system is fundamental to good maintenance. For 
this reason, the cameras are stored in a temperature- 
controlled office in their original cases when not in use. 
When being used during an experiment, these instruments 
are handled with the utmost of care. Of particular impor- 
tance is keeping the detectors cooled with liquid nitrogen. 
Before the cameras are turned on, the dewars are filled 
and the detectors are allowed adequate time to reach 
operating conditions. The dewars are continually topped- 
off with liquid nitrogen so that their response is stable. 

When the Agema system was purchased by NASA Ames 
in 1990, it represented the state-of-the-art in infrared 
imaging technology. During the time of this experiment, 
the state-of-the-art has seen the introduction of staring 
Focal Plane Array (FPA) systems. An FPA system 
typically incorporates an array of detectors as opposed to 
one single detector element (ref. 58). For this reason FPA 
systems provide image quality that far surpasses that of 
older-generation scanning systems. However, FPA 
systems currently suffer from calibration problems due 
to the large number of detector elements (ref. 59). 

The measurement accuracy of the Agema 880 infrared 
imaging system is maintained by constantly monitoring 
internal temperature sensors that compensate for detector 
response by adjusting the gain of the system after each 
scanning sequence. By “recalibrating” the detector at the 
end of each scan line, the older generation Agema 880 
system offers accuracy and predictability that cannot be 
attained by current generation FPA systems. Because the 
Agema 880 system offers this advantage over FPA 
systems, and because of its immediate accessibility to our 
testing crew, it was the logical choice for instrumentation. 

3.12 Uncertainty of Measured Infrared Images 

As mentioned in Section 3.11, the Agema 880 Dual Band 
scanning infrared camera system has a source sensitivity 
of 0. 1 °C at 30.0 °C for the MWB camera. The LWB 
camera has a source sensitivity of 0.05 °C at 30.0 °C. This 
uncertainty in the source sensitivity can produce errors in 
the measured radiance. Figure 49 shows how three color 
legends can be attached to the same measured infrared 
image by virtue of the instrument uncertainty. The three 
color legends show the minimum radiance due to 
uncertainty, the measured radiance, and the maximum 
radiance due to uncertainty, respectively. The scale for 
the minimum radiance is obtained by subtracting from 
the measured radiance the instrument uncertainty, and 


linearly scaling the color legend. The scale for the 
maximum radiance is similarly obtained by adding the 
instrument uncertainty to the measured radiance. In this 
example, the difference between the minimum and 
maximum measured radiance due to instrument 
uncertainty varies less than 0.27 percent. 

3.13 Midac Spectrometer Data 

One goal of the APU experiment was to gain a quanti- 
tative understanding of the levels of carbon dioxide, water 
vapor, and carbon monoxide in the APU exhaust plume. 
Charles Bennett of Lawrence Livermore National 
Laboratory (LLNL), in Livermore, California, offered to 
investigate with us the feasibility of using an infrared 
spectrometer for this purpose. Before the preliminary 
experiment we were hopeful that Bennett's technique of 
using the measured spectra to predict exhaust plume 
constituent levels would work. Unfortunately, this 
technology is not yet sufficiently mature, and so the focus 
of Bennett's post-processing analysis shifted to using the 
APU spectra to predict a pre-measured temperature within 
the exhaust plume. 

Bennett visited both the preliminary experiment in May 
1 993 and the dedicated experiment on July 1 3, 1 993, and 
took several hours of data with the Midac spectrometer, 
as described in Section 2.3.16. Post-experiment analysis 
of the spectrometer data yielded results that were 
unsuccessful in predicting the exhaust plume constituent 
distributions. The inability to resolve carbon dioxide and 
water vapor levels from the measured spectra is attributed 
to the fact that the spectra of these two species overlap 
significantly in the MWB. Because the spectra are 
inseparable in the MWB, it may be impossible to predict 
the levels of carbon dioxide and water vapor in the plume. 

Because carbon monoxide emits in a region of the MWB 
that does not overlap with the spectra of water vapor and 
carbon dioxide, the focus of Bennett's post-processing 
efforts changed. It was now the goal of his research to use 
the carbon monoxide spectra to predict the temperature of 
a particular point in the APU exhaust plume. By matching 
the measured spectrum with the spectrum corresponding 
to a known temperature, the temperature of the source 
could be determined. The digitally controlled traverse rig 
was used to position the probe along the linc-of-sight 
connecting the spectrometer with the light source. The 
temperature and pressure at that point were then measured 
by the probe. Each spectrometer survey was conducted at 
one diameter downstream of the nozzle exit. 

Predicting the temperature of the exhaust plume using the 
measured spectra is complicated by virtue of looking 
through a column of gas through which there exists a 
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steep temperature gradient. In this scenario, the 
spectrometer looks through the relatively cold ambient 
air, into the hot exhaust plume, back into the relatively 
cold ambient air, and finally into a light source having a 
known spectrum. An added benefit of conducting the 
surveys near the nozzle exit, one diameter downstream, is 
that the temperature profile in the exhaust plume at these 
locations is fairly flat so that the measured temperature 
would approximate the average temperature of the 
exhaust plume along the line-of-sight between the light 
source and the spectrometer. 

The process of predicting a temperature by matching the 
measured spectrum with a spectrum from a source of 
known temperature is illustrated in figure 50. Here the 
solid symbols represent the measured data and the open 
symbols represent the spectrum of carbon monoxide at 
some temperature. In figure 50, the spectrum correspond- 
ing to a known temperature approximates the measured 
data; when the two spectra match, the temperature of the 
emitting constituent is hypothesized to be that of the 
theoretical spectrum. The results from this analysis arc 
summarized in table 5. As can be seen from the results 
given in table 5, the temperature determined in this way is 
within ten to twenty percent of the measured temperature. 

4, Infrared Image Prediction 

Recall that this work is partially motivated by the need to 
validate computational fluid dynamic (CFD) codes. Also 
recall that the validation technique requires a quantitative 
comparison between an experimental infrared image and 
an infrared image predicted from a CFD solution. The 
latter image incorporates the use of an inverse Monte 
Carlo ray-trace technique. This section describes the ray- 
tracing technique. Infrared images predicted from the data 
obtained during the axisymmetric auxiliary power unit 
(APU) experiment are also shown. 

It is the intent of this report to provide to the scientific 
community a meticulously obtained and completely 
documented data base for infrared analysis of sub-sonic 
jet engine exhaust plumes. One motivation for creating 
the data base is to validate Turk’s (ref. 25) infrared image 
prediction algorithm. Improvements arc still being 
incorporated into Turk’s computer program, and so this 
report describes the current state of a work-in-progress. 
Turk’s doctoral dissertation describes the development of 
the Monte Carlo ray-trace technique. Specifically, his 
contribution focuses on methods for accelerating the ray- 
trace process. The reader is referred to Turk’s dissertation 
for a more detailed description of the Monte Carlo ray- 
trace method (ref. 25). 


4.1 History 

A virtual panacea of techniques is available for infrared 
image prediction of exhaust plumes. In general this 
technology has been developed for use by the defense 
community. As a consequence of limited availability due 
to security concerns, no single scientist or engineer is 
likely to have seen all of the methods that have evolved 
for predicting infrared images. However, over the past 
four years the scientific community has witnessed a 
dramatic improvement in the ability to predict infrared 
images that approach measured data. By comparing 
General Electric’s infrared image prediction program, 
LOIR (ref. 33), which was supposedly the state-of-the-art 
in 1990, to the images that can now be predicted using 
Turk’s inverse Monte Carlo ray-trace computer program, 
it is obvious that a significant contribution has been made 
to this field. It is now possible to attach a similar radiative 
flux scale to both the predicted and the measured infrared 
images, a capability nearly unimaginable only four years 
ago with LOIR. It should be mentioned that many 
improvements to LOIR have been made since 1990 and 
that the current version General Electric’s infrared 
analysis code is now called LOIR II (ref. 60). 

4.2 How the Ray -Trace Method Works 

The reverse Monte Carlo ray-trace computer program has 
been coded to predict the infrared image that would be 
captured by an infrared camera. The word “reverse,” or 
“inverse,” is attached to the name because the rays are not 
traced forward from the scene to the infrared camera as 
the phenomenon actually occurs. Rather, the rays arc 
traced in reverse from the camera to the scene. By tracing 
the rays from the observer to the scene in the reverse 
manner, only that radiation leaving the scene and arriving 
at the infrared camera need be modeled. Accounting 
for only this radiation significantly reduces the work 
required to solve the problem. As they are described in 
Section 4.3, the problem can be solved in the reverse 
fashion because of the reciprocity relationships which 
govern the radiative exchange of energy between surfaces 
and volumes. 

In order to develop a two-dimensional infrared image 
from a three-dimensional CFD or CFD-likc solution, 
several rays are traced from the camera detector clement, 
through a pixelized virtual screen, and finally into the 
object space or the solution domain. Figure 5 1 illustrates 
this methodology. The quality of the predicted image 
improves by increasing the number of pixels on the 
virtual screen and by increasing the number of rays used 
per pixel, within the resolution of the CFD solution itself . 
The predicted infrared image has converged when the 
object space has been adequately sampled such that an 
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increase in the number of rays does not improve the 
quality of the predicted image. 

When a ray leaves the detector element of the camera and 
traverses a pixel on the virtual screen, its direction is 
selected at random as shown in figure 52. Two random 
numbers are generated which define a point on the 
detector element. Next, two additional random numbers 
are generated which define a point within a given pixel 
located on the virtual screen. The orientation of the 
outgoing ray is defined along a line connecting the point 
on the detector element with the point within the pixel. 
Note that in this scenario both the detector and the pixel 
have finite areas. 

Once the ray leaves the infrared camera and pixel i zed 
virtual screen, it enters the CFD solution space. Here the 
ray encounters a cell, where a cell is a region of the CFD 
solution domain that is bounded by eight nodes. The ray 
can either be absorbed in that cell, be transmitted through 
it, or be scattered. Figure 53 illustrates all possible 
outcomes for a ray that enters a CFD solution domain. If 
the ray encounters a solid surface it can be absorbed by 
that surface or reflected from it. If the ray is transmitted 
through a cell, it continues to traverse the CFD solution 
domain until it is absorbed in another cell or until it exits 
the solution domain. The termination point of the ray 
determines the radiative contribution of that ray to a pixel 
on the virtual screen. 

When a ray enters a cell of the CFD solution several 
decisions must be made to determine if the ray is 
absorbed in that cell. To begin with, the eight bounding 
nodes of the cell are used to define the six surfaces which 
describe that cell, as shown in figure 54. The ray’s 
entrance point on one surface of the cell and its direction 
are used to compute its exit point. When the exit point has 
been defined, the ray’s path length, s , through the cell can 
be computed by 

J = “■*./')* +(>'/ ->'jY +(zi -ijf (6) 

where jc, y, and z represent Cartesian coordinates, and the 
subscripts i and j distinguish between the entrance and 
exit points, respectively, as shown in figure 54. 

Next, the average temperature, average pressure, and 
average species concentration of the cell are computed 
using Shepard’s interpolation scheme, which is described 
in reference 61 . Here, the average parameters are com- 
puted at the centroid of the cell by using the parameter 
values at the cell’s eight bounding nodes. The distance 
between the cell centroid and each bounding node is used 
to weight that node’s contribution to the final average. 


The average value, S(P^), of the temperature, pressure or 
species concentration at the cell centroid, P/, where P, is 
defined by is given by 

N=8 „ 


(7) 

where the summation occurs over the eight nodes, 
denoted by the subscript j , that define the cell, and the 
magnitude of the distance between the cell centroid P/, 
and some point Pj, where Pj is defined by (xj>\j,Zj) is 
given by 

h - P A = Jfc - x j ) 2 + (>«■ - yj f + ( z i ~ z jf (8) 

or 

( 9 ) 

which follows from equation 6. 

Once the average temperature, average pressure, species 
concentration, and path length have been defined for a 
cell, the absorptivity, a, of that cell can be computed 
using an infrared band model. The infrared band model 
used to compute cx\(X) is described in Nelson’s Master of 
Science thesis (ref. 22). The character A, in parentheses, 
denotes the wavelength dependence of this computation. 
All the equations necessary to make this computation are 
given in reference 22. 

Next, a random number is generated and compared with 
the absorptivity of the cell to determine if the ray is 
absorbed in that cell or if the ray is transmitted through 
the cell. If the random number is less than or equal to the 
absorptivity, i.e., 

random number <a\(X) (10) 

then the ray is absorbed in that cell. If the random number 
is greater than the absorptivity, i.e., 

random number > (A) (11) 

then the ray is transmitted through the cell. The sequence 
of events described above is followed until the ray is 
absorbed in a cell, or until the ray exits the solution 
domain. Upon exiling the solution domain, the ray can be 
absorbed by an infinite extension of the ambient atmo- 
sphere, a background simulation, or a blackbody at a 
temperature of absolute zero, for example. 
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4.3 Energy Computation 

The ray tracing technique is useful in computing the 
radiative flux or energy arriving at the observer location 
from the scene. This computation is more easily described 
after first presenting the reciprocity relationships for 
radiative exchange between surfaces and volumes. The 
equation used to compute the energy arriving at the 
observer location can then be more easily understood. 

The reciprocity relationships are developed here because 
in this analysis radiation is emitted from surface and 
volume elements (cells), and absorbed by a surface 
element (the detector element). 


and a volumc-to-volumc energy exchange. 

D Vj,Vj J e A.(7- c Vj • TV* • A/ ■•' v; (-4*)) dA* 

Aproj. V) 

(16) 

= D V y V, J £ X ( X. C V 1 , T v . . /Vj • -v Vj ( A * ) j dA * 

Aproj. Vj 

In equations 15 and 16, C is the concentration of species. 

P is the pressure, and s is the path length of the volume. 

A * is the projected area of the volume. 


The reciprocity relationship, in the absence of a change in 
the index of refraction, for a radiative exchange between 
surface i and surface j given by 

£>S, Sj EA..5, (*- T S, ) A S . - D Sj Si £\ Sj (x T Sj j A Sj ( 1 2) 


In equation 12, D is the distribution factor, £ is the 
hemispherical emissivity, and A is the area of the surface. 
The subscript A, denotes spectral dependence, and the 
subscript S identifies the variable as a surface property. 
The independent variables A and T denote wavelength and 
temperature dependence, respectively. The distribution 
factor, Dsi is the fraction of rays emitted by surface / 
that are absorbed by surface j due to direct radiation and 
all possible reflections, refractions, and scattering events. 
The distribution factor for the radiative exchange between 
surface i and surface j is estimated by 


°S, Sj 
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( 13 ) 


In this development, the spectral emittance, £is related to 
the absorption coefficient, or, by 

ex(X.C,r.P..s)= (17) 

The absorption coefficient, a, is computed by an infrared 
band model. 

By invoking the reciprocity relationships given in 
equations 12-16, the energy arriving at the observer 
location, o, from all surfaces and volumes is given by 

dE 0 (k) = Y D oi 7, )i\ h {X. 7 ■ )Ajdk 


N v 
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where Sj is the number of rays emitted by surface i 
and absorbed by surface j, and is the total number of 
rays emitted by surface /. The distribution iactor is also 
equal to the ratio of energy emitted by surface i and 
absorbed by surface j to the total energy emitted by 
surface /, that is, 


%, Sj = 
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A similar reciprocity relationship can be written for a 
surface-to-volume energy exchange, 

D s„Vj ex (x. r Si ) As, 

= 4 D Vj,Sj J tx(x. c v . , T v . , s Vj (4* j j dA* 

Aproj. Vj 


x 4K,'xb (x. Tj)j £\(....)dA* dX 

A J 

( 18 ) 

where i'b.X(T) is the Plank blackbody radiation distri- 
bution function. Equation 18 may be rewritten after 
canceling and collecting terms to give the energy arriving 
at the observer location as 

Hall 

dE n (X)=^n D ok £ X (X. T 0 ) i'\ {X. T k ) A„ dX (19) 
k=l 


4.4 Line-of-Sight Method 

One disadvantage of the Monte Carlo technique is the 
time required to obtain broad-band high quality images 
Because the ray-trace must be conducted over small 
wavelength intervals, AA, broad band images are lime 
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consuming to obtain. A broad-band image is obtained by 
integrating across the wavelength band; the broader the 
wavelength band, the longer it takes to predict an infrared 
image. 

The Line-of-Sight (LOS) method for predicting infrared 
images is faster than the Monte Carlo technique because 
fewer rays are required to sample the solution domain. 

The LOS method is like the Monte Carlo technique in 
that it utilizes a graphical ray-trace method. The sequence 
of tracing rays begins by generating a random initial 
direction for the ray as described in Section 4.2. When 
a ray leaves the pixelized screen and enters a cell, the 
path length through the cell is computed and the eight 
bounding nodes arc used to compute the average tempera- 
ture, pressure, and species concentration of the cell. An 
infrared band model is then invoked to compute the 
transmissivity, ZA/b °f the cell. The transmissivity of the 
ray is the product of the transmissivity of each cell 
traversed by the ray 

J 

^ AX, ray = Y[ x *m m 

i-0 

The ray continues to the next cell if the transmissivity of 
the ray is greater than some small value, say 0.001; i.e., 

TaJUov > 0-001 (21) 

In the LOS method, the ray continues to be traced through 
the solution domain until the total transmissivity of the 
ray is less than some small value, say 0.001; i.e., 

Kray *0.001 ( 22 ) 

The energy contributed by the ray to the computed 
infrared image is given by 

dE rayA K ) 


n 

= ^X,i-I ^\,b * \,i 7/ ) A () 

i=I 


(23) 
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where the sum is over the number of cells traversed by 
the ray, the emissivity and the ray transmissivity is 
defined in equation 2 1 . In equation 23 the differential 
solid angle subtended by the ray can be written 


clu> ray 


/ pixel 

yv r 2 


(24) 


In equation 24, dA is the projected differential area of a 
pixel to the observer, N is the total number of rays traced, 
and r is the distance between the observer and the pixel 
on the virtual screen. 


It is important to note that in the LOS method the 
contribution of each cell encountered by the ray is 
included in the equation used to compute the energy 
arriving at the observer location. The Monte Carlo 
technique is different from the LOS method in that only 
the contribution of the terminating cell is included in the 
Monte Carlo analysis, as indicated in equation 19. 
Because the contribution of each cell is included in the 
LOS analysis, fewer rays are required to adequately 
sample the solution domain. Consequently, the LOS 
method requires fewer rays to achieve convergence. 


4.5 Increasing the Number of Rays Per Pixel Improves 
the Quality of the Predicted Images 

The ultimate goal of the reduction of the three- 
dimensional CFD solution domain to a two-dimensional 
infrared image is a pixel-by-pixel comparison of the 
predicted image with measured experimental results. 

Such a comparison would not only serve as a validation 
technique for the CFD solution, but would also suggest, 
although not prove, that the solution domain was 
adequately sampled by the ray-trace method. It is intuitive 
that increasing the number of rays used to sample the 
solution domain would increase the probability of 
achieving a converged and accurate image, within the 
limits of the CFD solution. However, due to the long run 
times required to achieve such a solution the majority of 
the images presented in this section are less-than- 
convergcd solutions. 

Figure 55 shows a comparison of three predicted infrared 
images. Each of the images is predicted using the LOS 
method and is computed across the 4. 1 pm to 4.8 pm 
wavelength interval. The appropriate system response 
function, shown in figure 18(b), has been incorporated in 
the analysis. Each of the images has 50 x 1 00 pixels and 
the same arbitrary color scale. The difference between 
these images is the number of rays per pixel used in their 
generation. Figure 55(a) took two rays per pixel to create. 
Figure 55(b) and figure 55(c) look twenty and forty rays 
per pixel to create, respectively. The visually apparent 
difference between figure 55(a) and figure 55(b) is 
greater than the rather subtle visual difference between 
figure 55(b) and figure 55(c). The images appear 
smoother as the number of rays per pixel increases 
because the variability of the current pixel value from the 
converged pixel value decreases. These three images 
illustrate that the quality of the predicted image improves 
as the number of rays per pixel increases. 

The image shown in figure 55(a) took approximately 
two hours of CPU time to create using a Cray C-90 with 
128 megawords of main memory. The computation was 
completed on one of the eight processors, where each 
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processor has a peak operational speed of one gigaflop. 
The image shown in figure 55(b) is the result of averaging 
together ten images like the one shown in figure 55(a). 
Similarly to create the image shown in figure 55(c), 
twenty images like figure 55(a) were averaged together. 


4.6 Effect of Exhaust Plume Constituents on the 
Predicted Infrared Image 

As mentioned in Section 3. 1 ! the attempts to measure the 
exhaust plume constituents with the Midac spectrometer 
were unsuccessful. In an effort to use reasonable values 
for the exhaust plume constituents, concentrations for 
water vapor, carbon dioxide and carbon monoxide were 
extracted from reference 62. The concentration of water 
vapor at the core was assumed to be 0.03 percent by 
volume. The concentration of carbon dioxide at the core 
was assume to be 0.03 1 percent by volume. Finally, the 
concentration for carbon monoxide was assumed to be 
0.0 percent. These values are likely to be inaccurate to a 
certain degree. However, the predicted infrared images 
that are shown in subsequent figures show acceptable 
agreement with the experimental data. While we know 
that the concentrations of these constituents are inaccurate 
to some extent, we think their values arc at least close to 
the actual values based on the agreement between the 
predicted and measured data. 


4.6.1 Species Concentration Mapped to Total 
Temperatures- The core species concentrations, 
mentioned above, were mapped throughout the exhaust 
plume and into the surrounding atmosphere using 
Crocco’s similarity relationship (ref. 63). The analogy 
maps the spread of the core concentrations according to 
the total temperature. The concentration of a species at a 
point, Cp h is given by 


^ core Cgmb 

pi rp _rp 

1 n 1 n u 


where the subscript amb is an abbreviation for ambient 
and denotes ambient conditions, and the subscript o 
denotes a total value, i.e., total temperature. An assump- 
tion that the turbulent Lewis number of the exhaust plume 
is equal to one limits the applicability of equation 25. 


4.6.2 Carbon Dioxide Variance and Its effect on the 
Predicted Infrared Image- Figure 56 shows four infra- 
red images that were predicted using the LOS method. 
Each of the three infrared images have 50 x 100 pixels 
and were predicted at 4. 175 pm (2395 cm' 1 ). The wave- 
length was selected in a region of the spectrum where 
carbon dioxide is the predominant emitter. Figure 56(a) 
shows an infrared image predicted using the baseline 


constituents described above. Figure 56(b), figure 56(c), 
and figure 56(d) show images where the concentrations of 
water vapor and carbon dioxide in the core have been 
increased by ten, twenty, and thirty percent, respectively. 
Each of these images are shown with the same arbitrary 
color scale. The intensities of the predicted infrared 
images increase as the concentrations of water vapor and 
carbon dioxide increase. However, the visual difference 
between the baseline image, figure 56(a), and the image 
where the concentrations of carbon dioxide and water 
vapor were increased by thirty percent, figure 56(d), is 
barely perceptible. 

Because images can often be misleading, it is useful to 
look at a plot of intensity change versus change in carbon 
dioxide concentration. Figure 57 shows a plot that 
compares the change of baseline peak intensity as the 
baseline level of carbon dioxide is varied. This plot 
illustrates nicely the effect on the predicted infrared 
image obtained by varying the level of carbon dioxide 
in the core of the exhaust. The baseline intensity was 
computed for an observer looking at the core of the 
exhaust at 4.175 pm (2395 cm' 1 ) using the LOS method. 

4.7 The LOS Image Compared to Filtered 
Experimental Data 

Figure 58 shows a comparison of a measured filtered 
infrared image with an infrared image predicted using the 
LOS method. The measured infrared image has been 
optically filtered to include data only within the 4. 1 to 

4.8 pm wavelength interval. The predicted infrared image 
is computed across the same interval, where the infrared 
band model has been invoked to provide a 5 cm -1 
resolution across the interval. A similar color scale has 
been used on both images so that a direct comparison 
between the two images can be made. 

The primary difference between the experimental image, 
shown in Figure 58(a), and the predicted image, shown in 
figure 58(b), can be seen downstream of the nozzle exit. 
Here, the predicted infrared image is cooler than the 
experimental infrared image. The most likely explanation 
for this discrepancy is an improper accounting for the 
exhaust plume species concentrations in the core. Also, 
by mapping the concentrations according to the total 
temperature, an additional degree of uncertainty has been 
added to the computation. Even though these images 
appear slightly different downstream of the nozzle exit, 
the favorable comparison between these two images 
represents the current state-of-the-art in infrared image 
prediction technology. 
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4.8 Monte Carlo Ray-Trace Image Compared to 
Experimental Data 

Figure 59 shows a comparison between an infrared 
image predicted using the Monte Carlo ray-trace method 
and an experimental filtered-MWB infrared image. 

Figure 59(a) shows an infrared image predicted at 
4. 175 pm (2395 citH ). The predicted image has 
25 x 50 pixels and five thousand rays were traced per 
pixel. The predicted image took approximately seven 
days to run on a Silicon Graphics Indigo (SGI) 2 
Extreme® Computer that utilizes a 150 MHz IP2 
processor and sixty-four megabytes of main memory 
(ref. 64). The experimental infrared image shown in 
Figure 59(b) is a filtered-MWB image that has 
140 x 70 pixels and represents one Field of data. 

At first glance the predicted and measured images shown 
in Figure 59 look quite dissimilar. However, the predicted 
image shown in Figure 59(a) is not a converged image. 

If the Monte Carlo ray-trace code had been run suffi- 
ciently long to predict a converged solution the image 
would resemble the image predicted by the LOS analysis 
shown in figure 58. To this date a converged solution has 
not been obtained by the Monte Carlo code using this data 
set due to the excessive CPU time required for such an 
analysis. 

5. Conclusion and Recommendations 

The goal of the body of work presented in this report is to 
provide to the user community a data base for infrared 
analysis of jet engine exhaust plumes. To meet this goal, 
the temperature and pressure profiles of an APU exhaust 
plume were surveyed. Corresponding MWB, filtered- 
MWB and LWB infrared images were obtained. These 
data have been published in a variety of conferences 
(refs. 23 and 24), as well as this NASA TM. Information 
on how users may obtain electronic versions of these data 
are given in this section. 

5.1 The Goal 

Infrared imaging has been shown to be a valuable tool for 
heated flow field visualization. The ability of an infrared 
imaging system to resolve both large- and small-scale 
flow structures makes it a useful tool for validating CFD 
codes. The validation technique consists of the direct 
comparison between experimentally obtained infrared 
images and images predicted from the CFD solution. 
Predicting infrared images from CFD solutions requires 
reducing the three-dimensional CFD solution to a two- 
dimensional infrared image. A reverse Monte Carlo 
algorithm has been developed to accomplish this. The 


data base described in the preceding sections has served 
to help validate this algorithm. 

The temperature and pressure profiles obtained from 
the APU experiment were assembled into a CFD-likc 
solution. The reverse Monte Carlo computer program was 
utilized to predict an infrared image using the CFD-Iike 
solution as input. The predicted infrared images were 
compared with experimental data in Section 4. The 
favorable comparison between the two images lends 
credibility to the computer program. 

5.2 The Experiment 

The APU experiment described in Section 2 was 
successful in meeting the goals discussed out in 
Section 5. 1 . Temperature and pressure surveys of the 
APU exhaust plume were made from the exit plane to ten 
diameters downstream. Much effort went into ensuring 
the accuracy and repeatability of the data. This effort 
included calibration of instrumentation, a thorough check 
out of the data acquisition system during the preliminary 
experiment, and duplicate temperature and pressure 
surveys. 

The work presented in this report is unique and different 
from other temperature and pressure surveys of exhaust 
plumes in that a corresponding set of MWB, MWB- 
filtered, and LWB infrared data were obtained. An 
aluminized Mylar backdrop was utilized to create a 
relatively cold and uniform backdrop for the infrared data. 
The scanning infrared cameras were recalibrated prior to 
the APU experiment. The infrared cameras were treated 
with the utmost care during the experiment to ensure the 
integrity and repeatability of their response. The result of 
these efforts is a data base of infrared images that 
complement the temperature and pressure data base. 

5.3 Publication of Results 

Different portions of the research presented in this report 
have been published in a variety of locations. Nelson’s 
Master of Science thesis describes the selection of an 
infrared band model for predicting infrared images 
(ref. 22). The work described in that thesis was presented 
at an SPIE conference by the author (ref. 23). Birckelbaw 
and Nelson show in an AIAA paper the utility of a 
scanning infrared imaging system to resolve the flow Field 
about several V/STOL aircraft (ref. 19). A subsequent 
SPIE conference paper (ref. 24) and a paper given 
recently at a JANNAF conference at Lockheed Missiles 
and Space Company (ref. 60) describe the body of work 
shown in this report. The publication of this research in a 
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variety of forums and at different conferences subjects the 
work to a review by professionals that are familiar with 
this technology. The data described here have been well 
received by the infrared analysis community. 

5.4 Public Access to Data 

This report, originally published as a doctoral 
dissertation, is currently being published as a NASA 
Technical Memorandum (NASA TM) (ref. 65) to make 
its access more readily available to the public. Also, 
data collected during the rectangular nozzle survey (not 
presented here) will subsequently be published as a 
NASA TM (ref. 66). Electronic versions of the tempera- 
ture, pressure, and infrared data may be obtained by 
contacting: 

Larry D. Birckelbaw 
Mail Stop 327-7 
NASA Ames Research Center 
Moffett Field, CA 94035-1000 


5.5 Recommendations for Future Work 

The most important contribution that could be made to 
the body of work presented in this report would be to 
accurately measure the exhaust plume constituents issuing 
from the APU nozzle. Currently the levels of the exhaust 
plume constituents are estimated and the proximity to the 
actual values is hypothesized to be close enough when the 
predicted infrared images compare favorably to measured 
results. Quantitatively measuring the APU exhaust plume 
constituents would be an important contribution to this 
data base. 

Finally, as discussed in Section 3.11, the Agcma 880 Dual 
Band scanning infrared camera system was a state-of-the- 
art infrared imaging system in 1990. If an experiment 
similar to this one were to be conducted in the future, a 
more modern imaging system would likely yield higher- 
grade data than that which could be attained with the 
Agema 880 system. For this reason, it is recommended 
that for future tests, other options for infrared imaging 
such as focal plane array systems be investigated. 




Appendix A 


Temperature and Pressure Data 
















































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































nit 7 


nn • (com) 


i « d « 

|SE^E9| 




pooMon 


pocWon 



0 » v ) 


( In .) 


RBE2BBI 

000 

0000 

- 4.50 

0 006 

0.00 

0 16 

0 978 

-4 00 

0.006 

0 16 

0 31 

0075 

-3 75 

0.006 

031 

0 47 

0 061 

-344 

0.006 

047 

063 

0970 

320 

0.004 

063 

0 70 

0961 

-313 

0004 

0 70 

004 

0.056 

-297 

0 003 

004 

1 00 

0 052 

-201 

0.001 

1 00 

1 25 

0.043 

266 

0.011 

1.25 

141 

0.034 

-250 

0073 

1 41 

1 56 

0040 

-242 

0.137 

1.56 

1 72 

0 009 

-234 

0 225 

1 72 

1 80 

0077 

-227 

0.335 

1 68 

1 05 

0 063 

-219 

0.452 

105 

2.03 

0 927 

-211 

0.660 

203 

2.11 

0035 

-203 

0.667 

211 

2 19 

0 657 

-1 95 

0.736 

210 

2.27 

0469 

- 1.06 

0780 

227 

234 

0264 

-1 72 

0 052 

234 

242 

0157 

- 1.66 

0002 

242 

2.50 

0073 

- 1.41 

0936 

250 

2 66 

0002 

-1 25 

0962 

266 

1 

2.01 

- 0.003 

•1 00 

0052 

201 

207 

-0 003 

-004 

0063 

207 

313 

-0 002 

-0 70 

0059 

3 13 

3 20 

-0 001 

-0 63 

0970 

320 

3 59 

-0 001 

-047 

0 066 

3 50 



-031 

0 074 




-016 

0 067 




0077 

0.961 

0067 

0.067 

0.947 

0 956 

0932 

0017 

0020 

0.913 

0 032 

0 965 

0.969 

0 039 

0.666 

0.767 

0 597 

0411 

0251 

0 136 

0 062 

0000 

0002 

0004 

0005 


0005 


















































































































































































































































































































































































































































































































































































































































































































































































































































































































Tabic A 1 3. Y-axis pressure data measured two diameters downstream; see figure 27(a) 









































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































Table A26. Y-axis temperature data measured five diameters downstream; see figure 30(b) 
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Table A28. X-axis temperature data measured five 
diameters downstream; see figure 30(d) 
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Table A32. X-axis temperature data measured six 
diameters downstream; sec figure 31(d) 
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Appendix B 
Weather Data 
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Speed 
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Table B3. Weather data lor July 14, 1993 
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Speed 
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Direction 
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Table B5. Weather data lor July 16, 1993 


Time 

Wind 

Speed 
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57.3 

95.4 

14.686 

06:30 

1.55 

. 

18 

57.7 

92.7 

14.688 

07:00 

0.45 

83 

57.9 

93.2 

14.688 

07:45 

5.35 

96 

58.7 

81.7 

14.690 

Time- 

Averaged 

Values 

3.81 

34 

58.8 

90.5 

14.685 


n 






































































































Appendix C 


Infrared Image Data File List 




Tabled. Infrared data obtained July 12, 1993 


Data Filename 

File Size 

Description 

Time 

runOb.swb 

runOb.lwb 

1 0 frames 

Plume against 
Mylar backdrop 

3:40 

runOe.swb 

runOe.lwb 

5 frames 

Plume against 
Mylar backdrop 

4:13 

nmlb.swb 

runlb.lwb 

1 0 frames 

Plume against 
Mylar backdrop 

4:51 

runle.swb 

runle.lwb 

8 frames 

Plume against 
Mylar backdrop 

5:24 

run2b.swb 

run2b.lwb 

10 frames 

Plume against 
Mylar backdrop 

5:29 

run2e.swb 

run2e.lwb 

1 0 frames 

Plume against 
Mylar backdrop 

6:04 

run3b.swb 

run3b.lwb 

10 frames 

Plume against 
Mylar backdrop 

6:05 

run3m.swb 

ain3m.twb 

10 frames 

Probe shown 
in the flow 

6:38 

run3e.swb 

run3e.lwb 

7 frames 

Plume against 
Mylar backdrop 

6:48 




Table C2. Infrared data obtained July 13, 1993 


Data Filename 

File Size 

Description 

Time 

clear.swb 
clear. Iwb 

S frames 

Mylar backdrop 
without plume 

2:20 a m. 

spectl.swb 
spectl .Iwb 

6 frames 

Plume against 
Mylar backdrop 

2:06 a.m. 

spect2.swb 

spect2.lwb 

6 frames 

Plume against 
Mylar backdrop 

2:12 a m. 

spect3.swb 

spect3.lwb 

6 frames 

Plume against 
Mylar backdrop 

2:38 a m. 

clearO.swb 

clearO.lwb 

6 frames 

Mylar backdrop 
without plume 

2:56 a.m. 

spect4.swb 

spect4.twb 

6 frames 

Plume against 
Mylar backdrop 

3:12 a.m. 

clear3.swb 

clear3.lwb 

6 frames 

Mylar backdrop 
without plume 

3:32 a.m. 

spect5.swb 

spect5.Kvb 

6 frames 

Plume against 
Mylar backdrop 

3:40 a.m. 

clear4.swb 

clear4.lwb 

6 frames 

Mylar backdrop 
without plume 

4:01 arr. 

spect6.swb 

spect6.lwb 

6 frames 

Plume against 
Mylar backdrop 

4:11 a.m. 

spect7.swb 

spect7.lwb 

6 frames 

Plume against 
Mylar backdrop 

4:27 a.m. 

clears .swb 
clears .Iwb 

6 frames 

Mylar backdrop 
without plume 

4:34 a.m. 

spect8.swb 

spect8.lwb 

6 frames 

Plume against 
Mylar backdrop 

4:52 a.m. 

clear6.swb 

clear6.lwb 

6 frames 

Mylar backdrop 
without plume 

5:06 a.m. 

spect9.swb 

spect9.!wb 

6 frames 

Plume against 
Mylar backdrop 

5:10 a.m. 

clear7.swb 

clear7.lwb 

6 frames 

Mylar backdrop 
without plume 

5:17 a.m. 




Table C3. Infrared data obtained July 14, 1993 


Data Filename 

File Size 

Description 

Time 

clearO.swb 

clearO.lwb 

6 frames 

Mylar backdrop 
without plume 

12:45 a.m. 

run4.swb 

run4.lwb 

6 frames 

Plume against 
Mylar backdrop 

1:18 a.m. j 

run5.swb 

run5.lwb 

6 frames 

Plume against 
Mylar backdrop 

1 :47 a m. 

run6.swb 

run6.fwb 

6 frames 

Plume against 
Mylar backdrop 

2:17 a.m. 

run7.swb 

run7.lwb 

6 frames 

Plume against 
Mylar backdrop 

2:51 a.m. 

run7m.swb 

run7m.lwb 

6 frames 

Plume against 
Mylar backdrop 

3:25 a.m. 

ain8.swb 

run8.lwb 

6 frames 

Plume against 
Mylar backdrop 

3:27 a.m. 

run9.swb 

run9.lwb 

6 frames 

Plume against 
Mylar backdrop 

4:06 a.m. 

run9b.swb 

run9b.lwb 

6 frames 

Plume against 
Mylar backdrop 

4:19 a.m. 

runl O.swb 
runlO.lwb 

6 frames 

Plume against 
Mylar backdrop 

5:00 a.m. 

runll.swb 

runll.lwb 

6 frames 

Plume against 
Mylar backdrop 

5:31 a.m. 

njnl Ib.swb 
runl Ib.lwb 

6 frames 

Plume against 
Mylar backdrop 

5:48 a.m. 

runllm.swb 

runllm.lwb 

6 frames 

Plume and 
Survey Probe 

6:12 a.m. 

run12.swb 

run12.lwb 

6 frames 

Plume against 
Mylar backdrop 

6:14 a.m. 

run13.swb 

run13.lwb 

6 frames 

Plume against 
Mylar backdrop 

6:41 a.m. 
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Table C4. Infrared data obtained July 15, 1993 


Data Filename 

File Size 

Description 

Time 

clearO.swb 

clearO.lwb 

6 frames 

Mylar backdrop 
without plume 

1:06 a.m. 

run14.swb 

run14.lwb 

6 frames 

Plume against 
Mylar backdrop 

2:02 a.m. 

run15.swb 

run15.lwb 

6 frames 

Plume against 
Mylar backdrop 

2:35 a.m. 

run16.swb 

run16.lwb 

6 frames 

Plume against 
Mylar backdrop 

2:57 a.m. 

run22.swb 

run22.lwb 

5 frames 

Plume against 
Mylar backdrop 

3:28 a.m. 

run 22h.swb 
run22h.lwb 

6 frames 

Plume above 
Mylar backdrop 

3:29 a.m. 

run23.swb 

run23.lwb 

5 frames 

Plume against 
Mylar backdrop 

3:57 a.m. 

run23h.swb 

run23h.lwb 

6 frames 

Plume against 
Mylar backdrop 

3:58 a.m. 

run24.swb 

run24.twb 

6 frames 

Plume against 
Mylar backdrop 

4:26 a.m. 

run25.swb 

run25.lwb 

6 frames 

Plume against 
Mylar backdrop 

4.59 a.m. 

run25h.swb 

run25h.lwb 

6 frames 

Plume above 
Mylar backdrop 

4:59 a.m. 

run17.swb 

run17.lwb 

6 frames 

Plume against 
Mylar backdrop 

5:39 a.m. 

runlS.swb 

run18.lwb 

6 frames 

Plume against 
Mylar backdrop 

6:09 a.m. 

run19.swb 

run19.lwb 

6 frames 

Plume against 
Mylar backdrop 

6:44 a.m. 




Table C5. Infrared data obtained July 16, 1993 


Data Filename 

File Size 

Description 

Time 

clearO.swb 

clearO.twb 

6 frames 

Mylar backdrop 
without plume 

12:59 a m. 

noshield.swb 

noshield.lwb 

6 frames 

Plume with no 
foreground shield 

1:02 a m. 

shield. swb 
shield.lwb 

6 frames 

Plume with 
foreground shield 

1:31 am. 

run26.swb 

run26.lwb 

6 frames 

Plume against 
Mylar backdrop 

1:38 a m. 

am27.swb 

run27.lwb 

6 frames 

Plume against 
Mylar backdrop 

2:23 a.m. 

run28.swb 

run28.lwb 

6 frames 

Plume against 
Mylar backdrop 

3:29 a m. 

run29.swb 

run29.lwb 

6 frames 

Plume against 
Mylar backdrop 

4:02 a m. 

ain30.swb 

run30.lwb 

6 frames 

Plume against 
Mylar backdrop 

4:26 a m. 

run31.swb 

run31.lwb 

6 frames 

Plume against 
Mylar backdrop 

5:05 a.m. 

run32.swb 

run32.lwb 

6 frames 

Plume against 
Mylar backdrop 

5:50 a.m. 

run33.swb 

run33.lwb 

6 frames 

Plume against 
Mylar backdrop 

6:27 a.m. 

clear45.swb 

clear45.lwb 

6 frames 

Cameras and 
Mylar backdrop 
rotated 45° 
without plume 

7:25 a.m. 

set 45. swb 
set45.lwb 

6 frames 

Cameras and 
Mylar backdrop 
rotated 45° 
with plume 

7:29 a.m. 
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Table 1. Preliminary experiment data from May 14-15, 1993 


Traverse 

Data Plane 

Data Filename 

Comment 

Exit Plane Data 

y-axis 

x = 0.0, z = 0.15 

prelim, run 1 

negative y-axis 

y-axis 

x = 0.0, z = 0.15 

prelim, run 2 

repeat run 

j 0.5 Diameters Downstream 

y-axis 

x = 0.0, z = 2.5 

prelim, run 3 

negative y-axis 

y-axis 

x = 0.0, z = 2.5 

prelim, run 4 

repeat run 

y-axis 

x = 0.0, z = 2.5 

prelim, run 5 

repeat run 

y-axis 

x = 0.0, z = 2.5 

prelim, run 6 

repeat run 

y-axis 

x = 0.0, z = 2.5 

prelim, run 7 

repeat run 

| 1.0 Diameters Downstream 

y-axis 

x = 0.0, z = 4.5 

prelim, run 8 

negative y-axis 

y-axis 

y = 2.25, z = 4.5 

prelim, run 9 

positive x-axis 

| 2.0 Diameters Downstream j] 

I x-axis 

y = 2.25, z = 4.5 

prelim, run 10 

negative y-axis | 

{ 3.0 Diameters Downstream || 

y-axis 

x = 0.0, z = 13.5 

prelim, run 11 

negative y-axis jj 

| 4.0 Diameters Downstream || 

I y-axis 

x = 0.0, z = 18.0 

prelim, run 12 

negative y-axis [ 

| 5.0 Diameters Downstream || 

1 y-axis 

x = 0.0, z = 22.5 

prelim, run 13 

negative y-axis | 

| z-axis sweep to 2.2 Diameters | 

z-axis 

x = 0.0, y = 2.25 

prelim, run 14 



X 



Table 2. APU data measured from the exit plane through three diameters downstream 
(these data obtained July 12-16, 1993) 


Traverse 

Data Plane 

Data Filename 

Comment 

Exit Plane Data || 

y-axis 

x = 0.0, z = 0.15 j 

run 1 

negative y-axis 

y-axis 

x = 0.0, z = 0.15 

run 2 

repeat run 

x-axis 

y = 0.0, z = 0.15 

run 3 

complete x-axis 

x-axis 

y = 0.0, z = 0.15 

run 4 

positive x-axis 

| 0.5 Diameters Downstream || 

y-axis 

x = 0.0, z = 2.25 

run 5 

negative y-axis 

y-axis 

x = 0.0, z = 2.25 

run 6 

repeat run 

x-axis 

y = 0.0, z = 2.25 

run 7 

complete x-axis 

x-axis 

y = 0.0, z = 2.25 

run 8 

positive x-axis 

( 1.0 Diameters Downstream || 

y-axis 

x = 0.0, z = 4.5 

run 9 

negative y-axis 

y-axis 

x = 0.0, z = 4.5 

run 10 

repeat run 

x-axis 

x = 0.0, z = 4.5 

run 11 

complete x-axis 

x-axis 

x = 0.0, z = 4.5 

run 12 

positive x-axis 

| 2.0 Diameters Downstream || 

y-axis 

x = 0.0, z = 9.0 

run 13 

negative y-axis 

y-axis 

x = 0.0, z = 9.0 

run 14 

repeat run 

x-axis 

x = 0.0, z = 9.0 

run 15 

complete x-axis 

x-axis 

x = 0.0, z = 9.0 

run 16 

positive x-axis 

1 3.0 Diameters Downstream f 

y-axis 

x = 0.0, z = 13.5 

run 17 

negative y-axis 

y-axis 

x = 0.0, z = 13.5 

run 18 

repeat run 




Tabic 3. APU data measured from three diameters downstream through ten diameters downstream 
(these data obtained July 12-16, 1993) 


Traverse 

Data Plane 

Data Filename 

Comment 

3.0 Diameters Downstream (cont.) 

|| x-axis 

x = 0.0, z = 13.5 

run 19 

complete x-axis 

4.0 Diameters Downstream | 

y-axis 

x = 0.0, z = 18.0 

run 20 

negative y-axis 

y-axis 

x = 0.0, z = 18.0 

run 21 

repeat run 

x-axis 

x = 0.0, z = 18.0 

run 22 

complete 

1 5.0 Diameters Downstream 

y-axis 

x = 0.0, z = 22.5 

run 23 

negative y-axis 

y-axis 

x = 0.0, z = 22.5 

run 24 

repeat run 

x-axis 

x = 0.0, z = 22.5 

run 25 

complete x-axis 

1 6.0 Diameters Downstream 

y-axis 

x = 0.0, z = 27.0 

run 26 

negative y-axis 

y-axis 

x = 0.0, z = 27.0 

run 27 

repeat run 

x-axis 

x = 0.0, z = 27.0 

run 28 

complete x-axis 

j 7.0 Diameters Downstream 

y-axis 

x = 0.0, z = 31.5 

run 29 

negative y-axis 

1 8.0 Diameters Downstream 1 

y-axis 

x = 0.0, z = 36.0 

run 30 

negative y-axis | 

| 9.0 Diameters Downstream 

y-axis 

x = 0.0, z = 40.5 

run 31 

negative y-axis I 

!j 10.0 Diameters Downstream 

y-axis 

x = 0.0, z = 45.0 

run 32 

negative y-axis 
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Table 4. Results from calibration of the Kiel probe 


Pressure 
maintained 
by the 
calibration 
unit 

Pressure 
measured 
by the 
kiei Probe 

(PSig) 

(psig) 

0.000 

-0.0004 

0.100 

0.0986 

0.199 

0.1986 

0.300 

0.2992 

0.399 

0.3985 

0.500 

0.5016 

0.599 

0.5994 

0.700 

0.7028 

0.800 

0.8036 

0.900 

0.9074 

1.000 

1.0034 

1.100 

1.1047 

1.200 

1.2049 


-1.40 


- 0.20 


-0.27 


-0.125 


0.320 


0.067 


0.400 


0.450 


0.822 


0.340 


0.427 


0.408 









































































Figure 1. (a) Photograph of the McDonnell Douglas Harrier YA V-8B, and (b) MWB infrared image of the Harrier in a 30-ft 
(9. 14-m) altitude and 30-kt (15.4-m/s) velocity flight condition. 
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Figure 3. (a) Photograph of the Boeing 747 Space Shuttle Carrier Aircraft, and (b) MWB infrared image of the Boeing 747. 
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Figure 4. Predicted infrared spectral image at 4.3 pm of the Boeing747 Space Shuttle Carrier Aircraft. 




Figure 5. Photograph of the experimental equipment taken during the APU test at the OARF. 
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Figure 6. The large two-dimensional traverse used during the APU experiment at the OARF. 
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Njote : all dimensions in in, (cm) 

(not to scale) 


Figure 10. The flow straightening assembly (ref. 20). 
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Figure 12. Combination thermocoupie and Kiel probe. 
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Figure 15. Flow chart of the Labviev v® data acquisition and analysis routine. 
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Figure 22. An example of conditioned temperature data showing “ clipping . 
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Figure 24. Experimental temperature and pressure measurements obtained at the exit plane, z = 0.0 in 
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Figure 26. Experimental temperature and pressure measurements obtained one diameter downstream, z = 4.5 in 
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Figure 27. Experimental temperature and pressure 
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Figure 28. Experimental temperature and pressure measurements obtained three diameters downstream, z = 13.5 in. 
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Figure 30. Experimental temperature and pressure measurements obtained five diameters downstream, z - 22.5 in. 
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(a) (b) 

Figure 32. Experimental temperature and pressure measurements obtained seven diameters downstream f z = 31.5 in. 



(a) (b) 


Figure 33. Experimental temperature and pressure measurements obtained eight diameters downstream, z = 36.0 in. 









(a) (b) 

Figure 34. Experimental temperature and pressure measurements obtained nine diameters downstream, z = 40.5 in. 




(a) (b) 


Figure 35. Experimental temperature and pressure measurements obtained ten diameters downstream, z = 45.0 in. 











(a) 



(b) 


Figure 38. (a) Velocity profile at three diameters downstream ; z= 13.5 in. (b) Velocity profile at four diameters 
downstream, z = 18.0 in. 
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(b) 

Figure 41. (a) Velocity profile at nine diameters downstream , z = 40.5 in. (b) Velocity profile at ten diameters downstream , 
z = 45.0 in. 
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Four "fields" of data, A, B, C, D, each consisting 
of 140 x 70 pixels, are "interlaced" together to create 
one "frame" of data of 140 x 280 pixels. 

Figure 42. Illustration of the construction of a “field” of data. 
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0.0 0.5 1.0 


Figure 43. Three consecutively measured fields of filtered MWB infrared data. 
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Figure 44. MWB infrared images of the APU exhaust plume. 
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Figure 45. Filtered MWB infrared images of the APU exhaust plume. 
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(a) (b) (c) 

Watts/sr/m/m 



11.45 26.26 41.07 


Figure 46. LWB infrared images of the APU exhaust plume. 
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Figure 47. Averaged MWB infrared images of the APU exhaust plume. 
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Watts/sr/m/m 



9.10 10.60 12.10 


(b) 


Figure 48. Infrared images of t he aluminized Mylar backdrop. 
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Maximum Radiance 
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Figure 49. The uncertainty of a measured infrared image. 
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Figure 50. Results from the Midac Spectrometer. 


133 




Detector 

(Observer) 


Figure 51. The three-dimensional CFD solution space is reduced to a two-dimensional infrared image by the ray-tracing 
process . 
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Figure 55. Illustration of the improvement in the quality of a predicted image as the number of rays per pixel increases. 
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<d> 


Figure 56. Predicted infrared images depicting the effect of varying the baseline carbon dioxide concentrations. 
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Figure 57. Effect on intensity for a given pixel by varying the baseline carbon dioxide concentrations. 


140 




Figure 58. Comparison between an infrared image predicted using the Line of Sight method and an experimental filtered - 
MWB image averaged over six frames of data. 
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Figure 59. Comparison between an infrared image predicted using the Monte Carlo ray-trace code and an experimental 
filtered-MWB field of data. 
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